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ABSTRACT 

An  algorithm  to  solve  systems  of  linear  equations  for  many  right  hand 
sides  Ax  =  b  where  A  is  a  symmetric  positive  definite  banded  matrix  is 
described  for  implementation  on  Illiac  IV.   Two  codes  are  discussed  —  one 
for  band  width  less  than  6h,    and  one  for  band  widths  ranging  between  6h 
and  127. 
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1.0   INTRODUCTION 

The  solution  of  systems  of  linear  equations  Ax  =  b  where  A  is  a 
symmetric  positive  definite  banded  matrix  occur  quite  commonly  in  aero- 
nautical engineering,  electrical  engineering,  photo gramme try,  and  other 
applications  that  give  use  to  matrices  exhibiting  a  large  degree  of  inter- 
element  dependency  over  some  banded  area  of  interest.   In  many  cases  the 
system  of  equations  is  required  to  be  solved  for  a  large  number  of  right 
hand  sides  which  usually  tends  to  be  a  function  of  the  computational  power 
of  the  computer  employed. 

The  parallel  computer  ILLIAC  IV  can  now  be  utilized  to  solve  such 
problems  with  a  marked  decrease  in  execution  time  per  case;  and  as  a  result 
will  permit  a  greater  number  of  solutions  to  be  generated. 


2.0  METHOD 

A  method  presented  in  [l]  is  utilized  in  this  algorithm  to  solve 

Ax  =  b  (1) 

where  A  is  an  n  x  n  positive  definite  band  matrix.   The  method  consists  of: 

1)  Decomposing  the  matrix  A  into  the  product  of  a  lower  triangular  matrix 
and  its  transpose. 

2)  Given  the  decomposition  of  A,  using  forward  and  backward  substitution 
routines  to  determine  the  solution  vector. 

2.1  DECOMPOSITION  OF  THE  BANDED  MATRIX 


If  A  is   a  positive   definite  matrix  of  band  width  m 


i.e.      a. .    =  0 


(  |i-j  I    >  m) 


then  there   exists    a  real  non-singular  lower  triangular  matrix  L  such  that 


LL     =  A  where   I       -  0  (i-j    >  m) 


(2) 


The  elements  of  L  may  be  determined  by  equating  both  sides  of  the  above 
equation,  as  in  the  standard  Cholesky  decomposition.   However  in  this  method 
it  becomes  possible  to  take  advantage  of  the  band  nature  of  the  matrix. 


L  is  given  by: 

i-1 
I..   =  (a.  .  -  I  A.,  A.,  )/£. . 


p  =  max  (0,  j-m) 
j  =  i+1,  . . . ,  i+m 


(3) 


i-1 

I.  .    =    (a.  .  -   J   £ 
n     ii   .  *• 

k=p 


ik 


2)2 


p  =  max  (0,  i-m) 


(h) 


2.2   FORWARD  AND  BACKWARD  SUBSTITUTION 


Substituting  (2)  in  (l)  the  following  is  obtained: 


LL  x  =  b 


t 


Now  let   L  x  =  y,    and  the   solution  of  the   set   of  equations  Ax  =  b   can  be 


determined  by  solving 

1 )  Ly  =  b   for  y  and  then 

2)  L  x  =  y   for  x. 

Taking  advantage  of  the  banded  nature  of  the   system,  the   elements   of  y  and 

x  are   found  as   follows : 

i-1  p  =  max   (0,    i-m) 

*i  ■  (bi "  J  £ik  V^ii  .    n 

k=p  i=0,    . . . ,    n-1 

p  =  min    (i+m,   n-1 ) 


P 

x.    =    (y,    -        I       K.    x.  )/*.. 

i  i        ,     .,,      ki     k        li 


k=i+l 


i   =  n-1,    . . . ,    0 


Once  the   decomposition   of  A  is    determined  as    discussed  in  the  previous 
section,   only  the   forward  and  backward  substitution  logic  need  be   employed 
to   solve  the   system  of  equations    for  various    right  hand  sides. 


3.0   IMPLEMENTATION 


3.1   STORAGE  SCHEME 


Consider  a  7  x  7  matrix  of  band  width  =  3: 


A  = 


r 


a00     a01     a02     a03 


a01      ail     \2     al3     \h 


a02   al2   a22   a23  a2U  a25 


a 


03     ^3     a23     a33     a3U     a35     a36 


a 


Hi\     \  £l\      r-  3.1 


11+     *2U     "3U       hk     aU5     "6 


a^_     a 


a„^      ar 


I 


25       35     %5       55       56 


a36    aU6     a56     a66 


The   storage  scheme   for  the   Illiac   IV  algorithm  requires  that   element    a. .  be 
stored   in  row   i,    PE  j-i.      That    is,   the  main  diagonal   is    stored  in  PE  0,   the 
first  superdiagonal   in  PE  1,    and  so   forth  where  the  uppermost   diagonal   is 
stored  in  PE  m  (m  <  6h) .      Naturally  since  the  matrix  is   symmetric,   only 
the  upper  triangle   of  A  is   stored  in  memory. 


PE 
0  12  3 

row  0: 


1: 

2: 

3: 

k: 
5: 
6: 


aoo 

aoi 

a02 

a03 

all 

ai2 

ai3 

alU 

a22 

a23 

a2k 

a25 

a33 

a3U 

a35 

a36 

%k 

%5 

%6 

X 

a55 

a56 

X 

X 

a66 


X  X 


CONVENTIONAL  STORAGE 


ILLIAC  IV  MEMORY 


Thus    for  band  width   <  61+,   the  A  matrix  requires   n  rows   of  PE  memory   and  m+1 
PE's,   where   again  n   is  the  order  of  the  matrix  and  m  the  band  width. 

When  the  band  width  exceeds  6k  then  two  rows  of  PE  memory  are  required 
per  row  of  A  (it  will  be  assumed  throughout  the  remainder  of  the  paper  that 
m  <   127).      The  storage   scheme   is   depicted  below  for  m  =  66,   n  =   500: 


Row  of  PEM\   0 
row  0: 


row  1: 

row  2: 
row  3: 

row  866: 
row  867: 

row  868: 
row  869: 

• 

row  872: 
row  873: 

row  87I+: 
row  875: 


•\ 


00 


a 


0,61+ 


T.1 


a. 


1,65 


a> 


a, 


a 


a 


a 


a. 


PE 


a 


01 


02 


a 


0,65 


12 


a 


1,66 


l0,66 


13 


ll,67 


a 


1+33,1+33         1+33,^31+         U33,U35 


a 


1+33,1+97       "1+33,1+98       "1+33,1+99 


a, 


]43l+,U3l+       "l+3l+,l+35       ^1+314,1+36 
1+314,1+98       al+3l+,l+99 


a, 


a, 


a 


1+37,1+37       "1+37,1+38       "1+37,^39 


ak99>l+99       x 


row  999: 


x 


1+36,1+36       "1+36,1+37       "1+36,1+38 

X  X  X  X 


62 


l0,62 


X 


a 


x 


a 


a 


a 


1+37,1+99 


63 


a 


0,63 


Row  0  of  A 


a. 


^,63  "1,61+ 


Row  1  of  A 


a, 


a, 


a. 


x 


X 


X 


X 


U33.U95       "U33.U961    ^  k33 

of  A 


h3k,k96       %3k,k9T\    Row  ,3, 

of  A 


U36.U98       %36,k99^    Row  ^ 

Of     A* 


} 
} 


Row  1+37 
of  A 


Row  1+99 
of  A 


Note  each  of  the  last  6U  rows  of  the  A  matrix  actually  only  require  a  single 
row  of  PE  memory  but  since  n  »  61+,  for  the  sake  of  programming  convenience, 
two   rows  of  PE  memory   are   allotted   for  each  row  of  A. 


The  b  vector  of  right  hand  sides    is    an  n-element   vector  which   is   stored 
across  the  PE's,    i.e.    in  example  1, 


0 


PE 
2  3 


aoo 

aoi 

a02 

X 

X 

X 

■n 

al2 

"l3 

X 

X 

X 

a22 

a23 

a2k 

X 

X 

X 

33 

a3U 

a35 

X 

X 

X 

%k 

%5 

X 

X 

X 

X 

a55 

X 

X 

X 

X 

X 

B:  b0       \       \       b3       b5       b5 


In   example  2     (n  =   500), 


B: 


PE 
2      62       63 


b0        bl        b2 
b6U    b65    h66 


b128  b129  b130 


•b62    b63 
'b126  b12T 
'b190  b191 


bi+U8  bl+^9  bli50"bU99*X         X 


The  b   vector  is   overwritten  by  the   solution  vector  x.      Similarly  the  matrix  A 

T 
is    replaced  by   its    decomposition   into  LL    ,    i.e.   by  the  L  matrix. 


3.2     ALGORITHM  OVERVIEW 


Based  on   an   input    flag  the   algorithm  will   either 
l)      perform  the   decomposition  of  A,    and  proceed  to  compute  the  solution  vector 
x  of  Ax  =  b,   or 


2)      assume  the   decomposition   of  A  is   input,    and  proceed  directly  to   compute 
the   solution  vector  x. 


A  is   already  decomposed 
Solve  Ax  =  "b   for  new 
rieht  hand  side. 


DECOMP 
Decompose 
A  into   LLT 


-  First  time 
routine  is 
called 


Solve   Ly  =  b 


T 

Solve  L  x  =  y 

for  x 


C 


-  Forward  substitution 


-  Backward  substitution 


RETURN 


J 


3.3  MATRIX  DECOMPOSITION 

The  matrix  A  is  decomposed  as  described  in  section  2.1  and  replaced  row 
by  row  by  the  new  matrix  L.   Although  the  programming  logic  required  to  decom- 
pose A  is  much  the  same  for  any  row  of  A  as  seen  by  the  flowchart  on  page  13, 
for  discussion  purposes  we  will  consider  the  following  segmentation  of  A: 


rows  0  to  m-1 
(first  m  rows  ) 


rows  m  to  n-m-1 
(n-2m  rows ) 


rows  n-m  to  n-1 
(last  m  rows ) 


Rows  0  to  m-1  1 

i-1    2  p 
For  i  <  m  I.  .   =    (a.  .  -  /  £..  ) 

ii   ki0  ik 

i-1 

A. .  ■  (a.  .  -   J   fc...  4.,  )/£.  .      j  =  i+1,  ....  i+m 
ji     ij   i=j_m  Jk  lk   1X 

The  first  row  (row  0)  is  a  special  case  and  is  computed  outside  the  main  loop, 
The  first  column  of  L  is  given  by 

^00  =     l/aoo 

£.^  =     a.^/&^  i   =  1,2,.  .  .  ,m 

iO  iO      00  '    '         ' 

T 

and  replaces  the   first   row  of     A.      Note  that   L     is   actually  being  stored  in 

PE  memory  since   it    is  much  easier  to   compute  rows  than  columns   on  Illiac   IV. 


PEO 

PEL 

PE2 

PE3 

00 

\o 

£20 

1 
30 

all 

a21 

a31 

aUl 

a22 

a32 

al+2 

a52 

Rows  1  through  m  of  the   decomposed  matrix  L  are  formed  as   follows:      (let 
i   =   index  of  row  currently  being  determined). 

1.  Create  M0DEA  ,    a  Boolean  variable,   to  be  used  as   a  mode  pattern  with  bits 
0  to  m-i   set . 

2.  SUM  =   0.0 

3.  Perform  the   following  for  k  =  0,    ....    i-1 

a.  Turn   on  PE's   0  to  m. 

b.  Multiply  row  k  by  I        -   call  this   result  PROD 

c.  Route  PROD   i-k  left 

d.  Set   leftmost   off  bit    in  M0DEA  pattern 


e.  Turn  on  PE's  based  on  MODEA  pattern 

f.  Sum  =  SUM  +  PROD  (in  enabled  PE's) 
k.      Row  i  =  Row  i  -  SUM  (PE's  0  to  m  are  on) 
5.   Turn  on  PE  0  only 


6.  Compute  l/l..      =  1/  /Row  i   (only  PE  0  enabled) 

7.  Enable  PE's  1  to  m 

8.  Compute  lx.    =   Row.  x  (l/£   )      J  =  1,  ...,  m 

ji      l       ii 

9.  Zero  fill  row  i  for  j  >  m 

*   It  is  useful  to  create  the  reciprocal  of  the  diagonal  elements  of  L,  since 
multiplication  is  much  faster  than  division  and  this  element  will  be  used 
again  in  the  forward  and  backward  substitution  logic. 


Suppose  A  is  a  7  x  T  matrix  of  bandwidth  3  and  the  third  row  of  L  (i  =  2) 
is  being  formed: 


PEO     PE1     PE2     PE3 


*00 

*10 

£20 

So 

hi 

*21 

Si 

hi 

a22 

S32 

alt2 

a52 

a 


66 


The  table  on  the  following  page  depicts  the  contents  of  special  variables  of 
interest  at  various  steps  in  the  algorithm: 


10 


Step 


Contents   of 


MODE  A 


PE  0 


PE  1 


PE  2 


PE  3 


PE  h 


2b 


3e 


3f 
3b 


3e 


3f 


E,    El 

Row  k   (k=0) 

lk 
PROD 

MODE  A 

PROD 

SUM 

SUM 

E,    El 

Row  k    (k=0) 

lk 
PROD 

MODE  A 

PROD 

SUM 

SUM 


'00 


20 
£00X£20 


1 


'20 


'20 


'11 


'21 


£11£21 


21 


20 


10 


'20 


*10X*20 


'20 


'20 


'20 


£30*£20 

0 

0 

0 

£30,£20 

0 

p  p 

£20   +£21      |  £30£20+£31£21      ,  £Ul£21 


I 


ROW  i(i=2)        a     -   I   I 


1 


i 


'30 


'20 


£30*£20 


i 


-^^k     ;    a32-^^3k-V1%2-^1- 


°52 


20 


l81 

£31 

£Ui 

0 

£21 

£21 

£21 

£ 

hx 

£31£21 

£Ul£21 

0 

1 

1 

1 

0 

£31£21 

£Ul£21 

0 

0 

£30£20 

0 

0 

0 

21 


MODE 


0 
0 


Row  i 


£22j  a22"   ^    £2k 
k=0 


(Same   as   in   step  h  above) 
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Step  Contents   of 


8  Row  i(i-2) 

of  LT 


PE  0 


'22 


PE  1 


PE  2 


PE   3 


1  1 

a32"   £   V2k     %2\ ^aUk£2k 

I  £  52     22 


22 


22 


PE  k 


T 
Row  2   of  L  £ 


22 


'32 


%2 


'52 


Rows  m  to   n-m-1 

The  procedure  outlined  for  rows   0  to  m  is    followed  with  the   following  exceptions: 

Step  1  -  Mode  A  =  bit   0  only  set. 

Step   3  -  k  =  i-m,    ...    i-1. 


Rows   n-m  to  n-1 

The  procedure  as  modified  for  rows  m  to  n-m-1  is  followed  except  that  in  previous  rows 
m+1  elements  of  L  were  computed,  and  for  these  rows  less  than  m+1  elements  are  com- 
puted.  For  row  i  only  n-i  elements  are  computed,  and  therefore  at  step  3e  a 
mode  pattern  called  SVMODE  is  "AND-ed"  with  the  MODE  A  pattern  in  determining 
"SUM".   For  the  first  n-m  rows  the  boolean  variable  SVMODE  consists  of  bits 
0  to  m  set  since  m+1  elements  are  being  formed,  and  for  the  last  m  rows  n-i 


elements  are  computed  so  bits  0  to  n-i-1  are  set  in  SVMODE.   Let  n=T,  m=3: 
Row   PE  0   PE  1   PE  2   PE  3 


Contents  of  SVMODE  when  row  i  formed 
PE  0   PE  1   PE  2   PE  3 


0 
1 
2 
3 
h 

5 

6 


00 


11 


22 


'33 


'kk 


55 
'66 


£ 


10 


21 


20 


31 


£ 


30 


\l 


32 

U2 

52 

%3 

£53 

£63 

'5U 

z6k 

0 

65 

0 

0 

0 

0 

0 

1 
1 
1 
1 
1 
1 
1 


1 
1 
1 
1 
1 
1 

0 


1 
1 
1 
1 
1 

0 
0 


1 
1 
1 
1 

0 

0 
0 


].. 


m 


jm  <  i  s   n-m-] 


'n-m  <'   i  <  n-1 
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In  the  first  n-m  or  h   rows,  m+1  or  k  elements  are  being  determined  and  bits 

0  to  m=3  of  SVMODE  are  set.   For  the  last  k   rows  T-i  elements  are  computed 

in  row  i,  e.g.  for  i=5  ,  two  elements  are  computed  and  bits  0  to  1  are  set  in 
SVMODE. 


13 


c 


DECOMPOSITION 


5 


Initialize 
i  =  0 
SVMODE  =  PE's 
0-m 


[NEWROW]       p.    Ik 


m  =  Bandwidth  of  Matrix 
n  =  Order  of  Matrix 


(  NXTROW  i 


Set  up  proper 
loop  and  mode 
controls. 


Loop 
i  =  0,  n-1 


MODEA  =  PE  0 


p  =  l-m 


Yes 


rows  m  to  n-m-1  


No 


SVMODE  = 
PE's  0  to  n-i 


No 


p  =  0 

MODEA  =  PE'j 
0  to  m-i 


-last  m  rows 
turn  off 
rightmost 
enabled  PE. 


-I first  m  rows 

begin  k  loop 
at  row  0 . 


p.  Ik 


Ik 


Loop 

k=p  to  i-1 


1.  PE's   0-m  enabled. 

2.  PROD  =  Row  kx    L 

3.  Route  PROD  left    i-k  PE's. 

k.    Set   leftmost   off  bit   in 
MODEA. 

5.  Enable  PE's   determined  by 
MODEA  A   SVMODE 

6.  SUM  =   SUM  +  PROD. 


k   <   i-1 


fc 


=  i-1 


Row  i  = 
Row  i   -  SUM 


INEWROW   )- 


1.  Enable  PE  0  only. 

2.  Compute   &..    =    >WT  (in 

PE  0  only). 

3.  Row  i  =  1.0/Row  i  =  l/£, 


li 


k.    Enable   PE's  1  to  m 

5.    Compute   a .  .    =  Row  i   x  lA.  . 

Zero, 
to   bk 


6.    Zero,  fill   row  i    (PE's  m+1 


NXTROW 


i  <  n-1 


p.  13 


C 


i  =  n-1 


SOLVEIT 


> 


P.  15 


see  note  p.  9 


Determine   solution  vector  by 
forward  and  backward  substitution. 
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3.U   FORWARD  SUBSTITUTION 

T    T 
Now  that  A  is  decomposed  into  LL  ,  LL  x  =  b.   Ly  =  b  must  be  solved 

for  y  by  forward  substitution  and  L  x  =  y  by  backward  substitution.   Since 

T 
L  is  actually  stored  in  PE  memory  at  this  point,  accessing  elements  of  L 

must  be  done  in  a  diagonal  fashion: 


'00 


10 


11 


'00 


'10 


11 


20 


'21 


'31 


20 


21 


'22 


'31    32 


%2 


'33 


'k3         hk 


'22 


'32   V 
'33   £U3 


'kk 


Conventional  storage  of  L 


Conventional  storage  of  LJ 


Row  0  of  L 


Row  1  of  L" 


Row  2  of  L" 


Row  3  of  L" 


Row  k   of  L" 


T 


Storage  of  L  in  ILLIAC  IV 
*  See  footnote,  page  9. 
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The  elements  of  y  in  Ly  =  b  are  given  by 

y.  =  (b.  -Y  i.vyj/t..      p  =  max  (0'  i-m) 
1   k=p  lk  k   "      1-0,  ...  n-l 

The  summation  in  the  formula  is  obtained  simply  by  summing  the  elements  formed 

when  the  i   row  of  L  is  multiplied  by  a  row  in  PE  memory  containing  the  last 

m  y.'s  calculated.   PE.  of  this  row  contains  y.  .  (j  =  1,  ...m)  except  for  PE^ 
i  J  i-J  0 

which  is  zero: 

PE  0      PE  1       PE  2    ....   PE  m 

Last  m  y  's        0       y         y      .  .  .  .   y 

i  l-l        x—d  l-m 

Row  i  of  L         %..  £.  .  .      £.._....£.. 

n      i»i-l      i»i-2  i,i-m 

Product  0       I.    .    _y.  _  I.    .    0y.  _...■£..   y. 

i,i-l  i-l   i,i-2  i-2       i,i-m,  l-m 


i-1 

The  sum  of  the  terms  in  the  final  row  yields    )  £..  yn  .   Note  that  for  i  <  m 

,  f     lk  k 
k=i-m 

some  of  the  values  in  the  first  row  would  simply  be  zero  (since  less  than  my., 

i-l  x  S 

have  been  computed)  and  the  sum  )      I.,    y.  is  obtained  in  the  same  manner. 

k=o  lk  k 

Finally   it    should  be  noted  that    i=0   is   a  special   case   and   computed  outside  the 
main  loop: 

y      is  moved  to  PE  1   in  the   row  containing  the  last  m  y.  ' s   and  the   algorithm 
proceeds  to   determine  y    ,   y_   up  to  y        .      The   elements   of  b   are  overwritten  by 
the  y  vector. 


3.5     BACKWARD  SUBSTITUTION 

Now  L  x  =  y  is  to  be   solved  for  x  where 

k=i+l  l   =  n-l,   n-2,    ...    0 


IT 


A  similar  procedure  is  employed  as  that  used  in  computing  the  y  vector  except 

T 
that  it  becomes  easier  to  access  elements  of  L  in  this  case  than  it  was  to 

obtain  elements  of  L  in  the  forward  substitution  phase  of  the  algorithm. 

T 
Again  the  last  m  x.  's  are  multiplied  by  row  i  of  L  and  summed  to  obtain 

the  summation  in  (  5  )  above. 

For  i=n-l 


x.  =  x  _  =  y  . /£  _    n 
1    n-1    n-1  n-1, n-1 


x  n  is  moved  to  PE  1  in  the  row  containing  the  last  m  x. 's  and  the  algorithm 
n-1  l 

continues  as  before  to  compute  x   ,  ...  up  to  x  .   The  x  vector  replaces  the 
y  vector  in  PE  memory. 


18 


f      SOLVEIT       J 


Forward  substitution 
Solve  Ly  =  b   for  y 


A 

(forloopL 


1.  Enable  PE  0 

2.  SVMODE  =  PE  0  only 

3.  Compute  yQ  =  "^0^00 
h.  Overwrite  b  with  y 
3.  Route  yp  right  1->RGR 


Loop 
i  =  1,  n-1 


Yes 


No 


1.  SVMODE  =  SVM0DE+  PE  i 

2.  Enable  PE's  0  to  i 

3.  TEMP  =   Contents  of 
RGR   (last    i  Y[i]'s) 

k.        RGX    =     i     -    VV.    mimhpy 


1.  Pick  up  row  i  of  L 

2.  Compute  L[i,k]*Y[k]. 


With  PE  i   enabled  store  routed 
Y[i]'s 

PE 
0       1  2  3        ...    I 


TEMP 
RGX 


0     yi-l     yi-2yi-3"-   y0 
i      i-1        i-2        i-3. . .    0 


use  register  RGX  to  diagonally 

access  rows  of  L. 


0 


XC 


k  =  P»  P+l»  • • •  i-l 
p  =  max  (0,  i-m) 


p.  19 


o 
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ROWSUM 

'Compute   SUM)y 

the   sum  of 
above 

roducts, 


-Sm  =  J     L.kyk 
k=p 


1.  Enable  PE  0  only 

2.  Compute  y .=(b . -SUM)*l/  .. 

l   l         n 

3.  Overwrite  b-  with  y^ 

h.    Enable  PE's  1  to  63 
5.  Bring  TEMP  to  RGA  and 

route  TEMP  right  1  ■»  RCP  I 


Yes 


No 


TEMP  =  Contents  of  RGR 
(last  m  y[i] 's ) 

RGX  =  RGX  +  1 


FORLOOP 


p.  18 


i  <  n-1 


PE 


RGR: 


0 

1 

2      ...M     M+1...63 

0 

yi 

y        ...y          0...    0 
1-1          1 -m 

i  =  n-1 


p.  20 


Forward  substitution  completed 
Ly  =  b  solved  for  y 
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Back  Substitution 

T 
Solve  L  x  =  y  for  x 


0 


BKLOOP 


1.  Enable  PE  0 

2.  SVMODE  =  PE  0  only 

3.  Compute  x     ,    -  y     ,  / £     -, 

*  n-1         n-1     n-1, 

n-1 

k.    Overwrite  y  _  with  x  _ 
n-1       n-1 
5.  Route  x„  j  right  1  ■»  RGR 


Loop 
j  =  l,n-l 


Yes 


No 


1.  SVMODE  =  SVMODE  +  PE  J 

2.  Enable  PE's  0  to  j 

3.  TEMP  =  Contents   of  RGR 

(last    j   x[j ]'s) 


With  PE  j    enabled  store   routed  xCj^s 
PE 
0  1  2        J 

TEMP  0        x  x        _,_  . . .    x     . 

I  n-j        n-j+1  n-1 


1.  Pick  up  row   i  =  n-l-j 
of  LT 

2.  Compute  LT[i,k]*x[k] 


J  I     o 

1 

2 

•   •    • 

n-1 

i    |    n-1 

n-2 

n-3 

•    •   • 

0 

e 


i+l,    i+2,    ...    p 
min    (i+m,  n-l) 


p.   21 
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BKLOOP 


p.    20 


0 


</      ROWSUM 
Compute  SUM^ 
the   sum  of 
above 
^products 


1.  Enable  PE  0   only 

2.  Compute  x.=(y .-SUM)*l/£. . 

11  n 

3.  Overwrite  v.    with  x. 

l  l 

k.    Enable  PE's   1  to   63 

5.    Bring  TEMP  to  RGA  and 
route  TEMP  right   1  ->  RGR 


.lejL 


No 


TEMP  =  Contents  of  RGR 
(last  m  X[i]'s) 


J_<  n-1 


C 


3   =  n-1 


RETURN 


3 


P 


SUM  =  .  l.L±*\ 


k=i+l 


RGR 


0 

1 

PE 
2      .  .. 

M        M+l    . 

.   63 

0 

x. 

l 

Vr- 

x.        0      . 
l+m 

..      0 

Back  Substitution   Completed 

T 
L  x  =  y  solved  for  x 
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1+.0   ALGORITHM  USAGE 

As  the  algorithm  attains  maximum  efficiency  when  m,  the  bandwidth  of 
the  matrix  is  less  than  6U,  two  ASK  codes  were  written: 

1.  BS0LVE6U  -  to  be  used  when  m  <  6k 

2.  BS0LVE128  -  to  be  used  when  6k  <   m  <  127 

The  call  to  these  procedures  should  take  the  following  form: 

CALL  BSOLVE  q*}   (<the  A  matrix>, 

<the  {  }  vector>, 

<the  order  of  the  system  (n)> 

<the  bandwidth  of  the  A  Matrix  (m)> 

<decomposition  flag>); 

All  parameters  are  passed  as  addresses  whose  contents  are  described  as  follows : 
1.   <the  A  matrix>  -  occupies  n  rows  of  PE  memory  with  each  row 

containing  the  upper  diagonal  of  the  symmetric  band  matrix.   See  section  3.1 
for  an  example. 

If  the  user  desires  to  display  the  matrix  A  or  retain  its  con- 
tents he  is  advised  to  do  so  before  the  call  to  BSOLVE  as  it  is  overwritten  by 

T 
the  L  matrix  (A  is  decomposed  into  LL  ). 

2.  <  the  b  vector  >  -  the  n  element  right  hand  side  vector  b  of  Ax  =  b   is 
input  to  the  algorithm  and  stored  across  PE  memory  as  described  in  section  3.0. 
Again,  the  contents  of  the  b  vector  are  overwritten  by  the  solution  vector. 

3.  <  the  x  vector  >  -  the  n  element  solution  vector  x  of  Ax  =  b  is  output  from 

the  algorithm  and  stored  in  the  same  manner  as  b  above. 

k.    <   the  order  of  the  system  >  -  n,  the  order  of  the  system  is  core  constrained, 
vhere  m  is  the  bandwidth  defined  below. 

Restrictions  on  n: 

BS0LVE6^:   2m+l  <.   n  <  (2000  -  user  core) 

BS0LVE128:        n  <  (1000  -  user  core) 
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5.    <  the  bandwidth  of  the  matrix  >      -  m,   the  bandwidth  of  the  matrix  is   the 
number  of  non   zero   elements    in   given  row  of  the  upper  triangle  of  A,    excluding 


the   diagonal.      For   example,   n=5,   m=2. 

m=2 


a00      ^01      a02 


all      a12      al3 


a22      a23     a2i+ 


a33     a3i+ 


1 


Vn-5 


6. 


L  %h  J         J 

<  decomposition   flag   >   - 

=  0     the   decomposition  of  A  will  be  performed  and  the   solution  vector  x 

will  be   determined  for  this   right  hand  side. 
=  1     the   decomposition   of  A  is    input    (has    already  been  calculated)    and 

the   algorithm  will  proceed  directly  to   evaluate  Ax  =  b    for   a  new 

right  hand  side. 


STORAGE 

STORAGE 

RELOCATABLE 

CONTENTS 

PARAMETER 

REQUIRED 

MODE 

ADDRESS 

DESTROYED 

<the  A  matrix> 

N  rows 

Skewed 

Yes 

Yes 

J.X.     b 

<the  v   vector> 

(n-1  mod 
+1  rows 

6k) 

Straight 

Yes 

v   v  b  vec 

Yes}  .    , 

mpu1 

<the   order  of  system>  1  word 
<the  band  width>  1  word 

decomposition   flag>      1  word 


Yes 
Yes 
Yes 


No    } 

No 
No 
No 


x  vector 
output 


In   accordance  with   ILLIAC   IV  subroutine   linkage   standards,   the   contents   of  ACARS 
0   and  1   as  well   as    $D32   -   $D63  are  preserved;   however  the  user  is   advised  that 
$D0   -   $D31  will  be  overwritten  by  BSOLVE. 


2k 


5.0   TEST  RESULTS 

Example  1: 

An  ALGOL  procedure  to  enable  validation  of  the  ASK  routines  discussed 
in  this  paper  was  extracted  from  [l].   The  matrix  shown  below  of  band 
width  =  2  was  suggested  as  a  test  case  in  [l]  with  the  order  of  the 
matrix  =  UO: 


5 

-4 

1 

- 

4 

6 

-4       1 

1 

-4 

6  -4 

1 

1 

-4       6 

-4 

1 

6 

-4 

- 

1 

-4 

1 

1 

-4 
1 

6 
-4 

-4 
5 

with  b  the  vector  of  right  hand  sides  =  (l,  0,  0  ...  0). 

The  results  between  the  correct  solution  as  published  by  Martin  and 
Wilkinson  and  the  BS0LVE6^  routine  compare  quite  favorably  as  seen  in  the 
following  table  on  the  next  page. 
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Example  2 

A  =  {a. .}  where  a. .  =  10  +  i 


a. .  =  a. .  =  1  +  .01  i-j 
=  0 


i  =  1, . . .  n 


i-j  |  <  m   "^  1  -  1,  .  .  .  n-1 

1-j  I  £  m  J  j  =  i+1 ,  .  .  .  n 


Let  n  =  10,  m  =  k 


A  = 


10   1.01   1.02   1.03   1.0U 

11    1.01   1.02   1.03   1.0^ 


15 


1.01 

1.02 

1.03 

1.01+ 

16 

1.01 

1.02 

1.03 

IT 

1.01 

1.02 

18 

1.01 

19 


b  =  (l,  0,  0  . . .  0) 
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2  ,  dV?  \ur><i^b-)<<\f     0  0 
l,i«2857U23l7t    00 

BSOLVE6I1  (ssk  simulator) 


S  $  K  X  irtVE:IcR  =  00?-Ju^  r, 
riaST  LOCATION  EXECUTED! 

THE    SOLilT jON    tfEOTOq 


00V346    0:16 


l  ncATio-j!  k      <"(  i.nr...Ti  rvo 

rfLOCATION    +nnOr    oi:i)  C(Ln<"ATlON    +0000    02:*)  r(Lnr,iTlnN    +0000    03:8 

♦  ,l«57l42**i7i<t!J0  iT»*  0001    +  ,?*^7  \H?-*57  lila7*l*>*n()0l         4..3ia?«S7i42fi621l3»*000l  +  .2JS7  1  a  ?  B  5  7"  14?  692a*  0001 

♦  ,2ia?*S7ia?aM217»*nnfti    +.  1  1  <*?db7 1 4?3^'d  3«**oooi         0  o 


.*,**•   *******   4   *   »   *   .    *****.  *l**i******.«»* 


*************************** *************' ******* «**********«*•****< 
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Example  3» 

The  same  A  matrix  and  b  vector  as  tested  in  example  two  was  tested  with 
BS0LVE128  with  n  =  76,  m  =  66.     In  most  cases  n»m,  but  due  to  the  slow 
speed  of  the  simulator  it  was  not  possible  to  let  n  become  as  large  as 


desirable. 
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APPENDIX  A  -  BSOLVE6U  ASK  LISTING 


htr.iN  oooooioo???? 

chws    6a>  ooooopoo???? 

TILL      lpHl  PPPP0300???? 

DEFINE  PRlNTAs  P0000400???? 

STA  TFf.Kp;  Lll(3)  Bl  pT?MPP*3»Tt  NPpj  PlSpLAyP  SCM^UJ)  00000*00???? 

8 DO  I     *  ZERO  00000 ft nO???? 

S,,1  ;     %    FIXED  PHIM  ONE  On00070U???? 

sup i    t   fixed  poim  64  oodooroo???? 

$o«j     *  Anr.f-Ess  Or  MTy  00000900???? 

*;)5:     J  AnHKSS  Pr  RHS  00 0  01000???? 

ii)ftj      *  INPE*  UF  fN-M)TH  PPW  Or  MjX  OOOO 1 1 00???? 

a    I.E.  FT^ST  OF    LAST  M  ROw$  0000 1  ?00???? 

,K|      C°L'       V)7;     *>  USED  IN  CREATING  MODE  eATTFt?NS  PPC01300???? 

i   vrTF:  J  LEO  K  LEO  h  PPOOlflOO???? 

•  Mi  E^U  Snfts  *  BANT  PlDTM  (  INPUT  }  00001500???? 
,N,  F.kU  S)9;  ?  DIMENSION  VF  MATRIX  00001*00???? 
.Nit  E^U  S|)l0!  X  N-l  00001700???? 
»FLAG*   E^U       Sol])    *  =o  FIrSt  OftLL  yO  POUTl»'E  OoOOlflOO???? 

•  PECJ  E^U  &  r)  1 5 1  *  MOPE  F'ATTFFN  FOR  FFO  T^'LY  OOOC 1  90C  ???? 
.PFHOUE-E^U       SMf.|    *  MnDF.  F^TTFPN  EH»  FF-S  r*  TO  «  POOOPOOO???? 

•  NOTPtOJEWU  Sn17»  X  MOrE  PATTERN  FOR  PE-S  1  Tn  M  00002100???? 
.SVMODE  iE«U  $Tlp:  X  SAVE  h'GDE  OOCOPPOO???? 
.MooEA:  E«U       Sol  95    X  0O00P3P0???? 

•  SAVEOj  EWJ       Sr>Po»    *  SAVf  ACAPp  OPOPPAPO???? 

•  SAVE3:  Etu  S,?3:  *  SAVE  AC*P3  0p00?500???? 
,SVLOUP:E«U       Sn2flt    x  SAVE  LOOP  INDEX  (0  TO  *'-l)  00002ftOO???? 

•  RETURJ  El*U  S.-)25j  *  RFTURf  ADOPESS  0O00P7OO???? 
.USP^DEjEWU       %■)?(,;          h     USFP  PDDE  PATTFRN  PPOOPPOO???? 

•  USRC°:  E«u  S,27s  X  ijSEr  AcAp  0  00CC2QO0???? 
.USRCli  E«U       $r>2f:    *  USER  ACAP  1  00003000???? 

•  (0003 1PU???? 
TEMPI  BLK  1;  X  TrMpCRAPY  ROw  STORAGE  0PO03?pO???? 
TEMPP»  (HLK  I;  *  TF^'P  STORAGE  00003300???? 
PEN!    'l.ATA  0.l»2,3»e»5»6»7»P»9»l0t1l»l?»i3»l«»1^»lf »17»1B»  PPP034Q0???? 

l9>?'.f?i,2?>23»?4,?5*P6.?7,?p.,2?»3o.»,1»3?,33»34,  00003500???? 

?S»36,3  7»3H»"»c'»A'#^l»'i?»'ti3»^^»^,5.ii6»'7»^P,/i9»50»  0000  3  ftp  0???? 

5i»52.cS^.5<4.55»5f»57»5fl.59.ftr,61»6?»fc3l  PO0  0  370O???? 

•  PPCQ3P00???? 
X  SUtRpLllNE  RO*"Sm  CALC'LATES  T;<E  SUM  Or  THE  FLE"ENTS  0PCC39PP???? 
X  IN  RPA.  RGA  ^  T  Ll  CONTAIN  THE  pESl'LT  Or  THE  SUM  UPON  OOOP^PPO???? 
X  EXITING  THE  ROUTINE  00004100???? 
X  O0004pP0???? 
ROWSMM  FlLLt  0000*300???? 

S 1 L ( 3 1     ,SAVE3»   *  SAVF  RETURN  APOR  COCCAttOC???? 

LlL(?)     .SVmODFj  X  riPF'E^T  MOHF  TO  ApAR  ?  00004*00  ???? 

(;L;KPf(2.)l            ^  PFVF«SE  pATTFPN  000 04 ft 0 ''???? 

LOEEl      T.p?i      X  TURK  ON  pE~S  NOT  DESTPEn  FOR  SUM  OpO0«  70<???? 

rLRA;             %    Arp  cLEap  Their  PGA  nnooaanr???? 

sete     F.pp.-r;  *  turn  dn  all  pe-s  c;r.ooft9oii???? 

SETFl      f.MD.rJ  f     Ff;P-LOGSUM  00005000???? 

%  00005100???? 

t_lT(2}     1»ft.lj    X  LPPF  COUNT  po005?OO???? 

SL-1T(3)    1»        X  SHIFT  IMPEX  00005300???? 

RTl        ?,A,n(3>;  X  ROUTE  ANn  00005^00???? 

AtPN       <Rf       *  ADD  00005500???? 

rSHL(3i    11        X  rOURLE  SHIFT  nIST  OP005ftOO???? 

TXLTAF'(2)  ,'bl                X  PDNTiNUF  00005700???? 
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END 


LDL(3>     .SAVE3I   *  GET  RETURN  Apo* 
E*CHL(3)   SlCPl     %    AND  RETURN 

pF    SUPPnUTl\jE    PDwSM 


BSOLVEfr/ilLNTPViM       FTLLl 


SU(3> 
SFTCO) 
SU(  3) 

SU(0) 

sUdj 

SFTE 
CLC(0)> 

SfL(O) 
L IT(0) 
S  U  (  o  > 

ilKoi 

SU  (01 
I  IT(0) 
SlL(0^ 


THF 
—  FL 


rn»M 
call 


r»" 

•  »»«i»f.iMi» 


L  I-'  A  0  (  2  > 
rSnR(O) 

STL  CO! 

ALIT(?1 

LbADr?i 
SU(0, 

ALIT(2l 
L0AD(2l 
LUA0(P1 
SlL(O) 

CSUB(0T 
SI L(0) 
ALIT(?> 
LUAP(?l 
L  U  A  0  (  0  > 
SUtO) 

LPLCO) 

C  5>  U  R  (  0  ) 
SU(0) 


.RFTitRj  x  SAVE"  PETuRN  AppRFSS 

E;       I  GrT  US£P  M0Df  PATTpRN 

•  i  j  S  R  h  D  E  t  2  Ak'D  SAVE  TT 

•  jiSRCOi  *  SAVF  ACAp  0 

•uSPr1 *  *  SA^r  *CA*  1 

E.OP."E'  SpTfl   E«*NP«E» 


•zero* 
■  1 1 

•  nNEi 

■  <S4|. 

.  s  I  y  t  4  j 

siooooo^oooooroooooooooisj 

.pEO,    »  MHr£  FUP  1  pE 

flp  TwE  'ALL  TO  PSOL"E  IS 
BcPLVE  'A.  F.  hi,  H,  FLAG,  j   WHFPF 
IS  AN.  T'OirATpP  sO  FNTEPJN",  ThF  PPUTTwE 
T-,1  fIR«?T  TTmF  fUEcnMpnsF  A),  OTMERwTSE 
Tr  SOLVE  ax=F  fpr  v,Ew  phS 
nr  axsf  anp  win  p^  ovfp- 


j jMP  TO  LOGIC 

is  TmE  "ATRTX 

WRITTEN  as  TT  IS  DECOMPOSED  INTO  L  K    [ 

U  T^E  :ATR7X  f  UF  AX=F  AMp  WILL  ALSH 

Be  OVERWRITTEN  PY  X  IN  THE  APOVF  EQN 

IS  tuE  -»ROER  HF  THF  MTX  i_EO  ?000 

IS  TME  aA\'D  «IDTH  PF  THF  MJX  LEO  64 


SCO! 
.  ADPA  J 

1) 

*C0: 

•  AnRrJ 

1  I 
SCO  I 

SCO  J 

•  Hi 

.<->ne  J 

1| 

*co» 

»C0| 

.HI 

•  M  I 

.lastm; 


%  AOPR  pF    MjX 

*  SmIfT  C|j  ador  TO  gET  R°W 

*  PE  AUPR  PF  MTX 

*  Bump  |,IST  POINTER 
a  An or  nF  F  vfctop 

%  STORE  rU  AnnR  IN  App 

*  BljMp  ACAR?  FOR  NEvT  PARAMETER  AODR 
1  PICK  UP  PARAMETER  IN  fO 

I  VALUE  IN  AnpR 

%  DIMENSION  OF  MTX 


I!  N-l 

*  BllMp  ArAR?  FOR  NEXT  pAOAMETFR  AQDR 

?!  PTCk  ijp  PARAMETER  IN  $f*0 

*  value   in7ahdr 
%  band  moth 


*  N-M 

%   Index  of  (n-m>th  pow  o*"  mtx 


0000* 
00005 
00  0  0  6 
00006 
C0CC6 

00006 
00006 
0  0006 
00006 
00006 
00006 
0  00  0  6 
00007 
00007 
00007 
00007 
00007 
00007 
00  00  7 
00007 
00007 
00007 
0000ft 
00008 
00008 
00008 

ooooa 

00008 
00008 
00008 
0000  » 
0000« 
00009 
00009 
00009 
00009 
00009 
00009 
00009 
00009 
000  0  9 
00009 
00010 
00010 
00010 
00010 
00010 
00010 
00010 

00010 
00010 
00010 

ooon 

00011 
00011 

ooon 

00011 


AOO???? 
900???? 
000???? 
100???? 
20G???? 

300???? 
400???? 

500777? 

600???? 
700???? 
fiOO???? 
900???? 
000???? 
100???? 
?00??7? 
300???? 
400???? 
500???? 
600???? 
700???? 
800???? 
Q  0  0  7  ?  7  7 

000???? 
100???? 
20-')???? 
300???? 
H  0  0  7  7  7  ? 

soo???? 

600???? 
700???? 
A0<<???? 
900???? 
000???? 
100???? 
200???? 
300???? 
4007777 

SOO???? 
600???? 
700???? 
AOO7777 

900???? 
000???? 
100???? 
200???? 
300???? 
4  01)77?? 
500???? 
600???? 

700???? 
800???? 
900???? 
000???? 
100???? 
200???? 
300???? 
400???? 
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LDLCO)    .Ml 

LUA       PEN| 

stTC(n>  i  j 

CUMPC(O)  I 
H'FEl     *CO 
SIL(O) 


.  s  V  I-  P  f)  E  >  % 


BAND  WIDTH 

PE  INplClES 

I  BITS  ON  WHERE  PFN>BA»'D 

MpVF  T  BITS  TO  ACARO 

REVERSE  PATTERN 

Pr-s  0  10  •"  ARE  ON 

SAVE  IN  Pr^OOE 

also  in  svmQdE 


WTr>TH 


LIM.  ^0">         .  PEC! 

CUMPC(?)|  *  PE-S    l-t>3    ON 

C**NO(<M        .pEnnoEi*  pF.-S    1     T0    m 

SlL(O)  .nOTpEO:*  SAVE    IN    NDTPEO 


A  L  I  T  (  2  i 
LtAO(?) 
H'ADd  i 
STL(i) 


l; 


%  AnDR  (,E  FLAG 
t  COME  MS  m  0  w 
*  SAVE  IN  APR 


PARAMETER 
IM  AfAR  1 


ZtHTA(l)  »+l;     %     IT  ELA3=0  THEN  DECOMPOSE  MTX 
jOMp      SoL^lTi*  OTHERWISE  jUMP 


BEGIN  MAIM  Lnno  -  REDUCE  A  Tp  Lp»'ER  TRIANGULA0 

-  iCn  COrTAlKS  AOPRESS  Or  ROw  TO 
r£  RFPLACFt'  (A  IS  Hr5TR0yE0> 

-  SC1  IS  LOOP  COUNTER  =0,1. ...N-l 

LlT(l)  0.1,01  %  SET  UP  ACAR  1  -  CDRRSP.  TO  T 

CADO(  n  .  sii  j     ;.  Tn  LOOP  FROM  0 

CKf'TLM)  2a;      2  UP  TO  *-l  (RY  nNE) 

STL(1)  .SVLOOP'*  SAVE    EUR    USE     IN    EOP'RA"K    S"RST    LOGIC 


LLLCO ) 
Ll-A 
Ji  HP 

* 

NXtRO"' :LOL(0) 
C^OD( 0) 
Li  L(2) 
LbEEl 
Lt>L(2) 

SU(?l 

LLL(3) 

CSUB( 3) 
LLSSK1  ) 


•adpaI  *  first  *nw  is  spec  case 

0(0}?  *  LOAD  IT  IM  RGA 

NEwRpwi  %    Ji'MP  OVER  MULTIPLYING  LnGTp,  ETq. 

.ADrA>  %    ADJUST  ACARO  TP  CPNTAI*'  Th* 


$ ,-  1  I     * 
•pFmOUE  t  % 


%r2t 

•  ONE» 
.<> 
Sf.l| 

•  Ml 


ADDRESS  Of  THE 

pe-s  o  to  m 
Turn  cn  pe-S  o 

SET  Ksl.   WILL 

when  i<m 
i 

I-f  U  *C3 


1-TH  RP''  OE  MT* 

TO  M 

ONLY  CHANGE 


S.)6,  rK  A'3'*  *  JUMP  IP  I  <■  (M-W) 


for  the  last  m  po"s  c'eate 
specific  mode  patterns 

ll'lc?)   ,svi'.pde?x  load  current  "noe  pattrrn 

C5HL(2}   1|       *  RIGHTMOST  ON  Pe  N0W  nrr 

LLEE1     S„2j     %  TijRN  OFF  PE  N-I  (   M,    M-1  ...   1) 

*  CPRPSP  VAL  OF  Is(N-M,N-(M-l ) , ,.N*n 

STL(2>    .sVmpDE»*  SAVE  NEW  MODE 
JUMP      FiXac2j 


000115 
0001  If, 

0P01 17 

oooilft 

00P1 19 
000120 
000121 
000122 
000123 
oooi 2a 

00012S 
000126 

000127 
00012ft 
000129 
000130 
000131 

000i32 
000133 
00013A 
000135 
0001 3f 
000137 
00013ft 
000139 
00014O 

o  o  o  1 4  i 

00014? 
0  0  014  3 
0001 4A 
000145 
000146 
0  0  0  1  A  7 
0  0  0  1  A  ft 
0001 49 
000150 

000151 
0001 5? 
000153 
000154 

000155 
0  0  0 1  5  6 

000157 
00  015ft 
000159 
000160 
0001 61 
000162 
000163 
0  0  016  4 
000165 
000166 
000167 
OOOlf ft 
0  0  016  9 
0  0  017  0 
000171 


00???? 
PC???? 

00???? 
00???? 
00???? 

00???? 
00???? 

oo???? 
oo???? 

00???? 

00???? 

op???? 

00???? 

on???? 

00???? 
00???? 

oo???? 

oo???? 

00???? 
00???? 

oo???? 

00???? 
00???? 
00???? 
00???? 

on???? 

00???? 

oo???? 

00???? 
00???? 

oo???? 

00???? 
00???? 
00???? 
00???? 
00???? 

00???? 
00???? 
00???? 

on???? 

Oy???? 
on???? 

00???? 
00???? 
00???? 
Oo???? 

0  0???? 

oo???? 

00???? 

or???? 

00???? 
00???? 
00???? 

po???? 
00???? 
00???? 
00???? 
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CMC3I   USSFU)  %n*,rlX\Z*tl    JUMP  IT  M  LEO  I  t  I*  <*-*>  O00l7?00???? 

*  00017300???? 
t  FOR  FIRST  SI  ROaS  -  ADJUST  VALuE  OF  K  ANO  SET  -  -  00017400???? 
I  L"HFR  LIMIT  OF  LnriP  «0  (R0V  0)  00017*00???? 
I  00017600???? 

LL»L(3)  •  it  %  H  00017700???? 

CSUR(3)  Sell  *  M.I  00017ROC???? 

STL(3)  ,<;  »  K  ■  M  -  i  00017900???? 

CLC(3)j  %  START  AT  FIRST  ROw  ooOlBOOO???? 

%  oooifijoo???? 

FIXAC21 ILITC2)  0,1,01  %  SFT  UP  ACAR  2  -  CORRSP.  TO  j  00018200???? 

CADD(?)  sell  <  TO  LOOP  TO  00018300V??? 

CSnB(2)  .oNEl  *  1-1  (tCl-1)      (UPPER  LlM^T)  00018400???? 

C^CTLC?)  2ul  %  FPOM  MAXO«I-M)  (LOWER  LIMIT)  0001*500???? 

CADDC2)  $r3l  00018600???? 

X  00018700???? 

I    THE  FOLLOWING  i_OnP  IS  TRA'FR^ED  M  TIMES  00018800???? 

1    EXCEPT  *'HEM  I<M  -  THFN  OMLY  T  TIMES  00018900???? 

%  00019000???? 

STL(O)  .SAVEO;  ?  SAVE  ACAR  0  00019100???? 

CLRAI  00019200???? 

STA  TEMPt  %  ZFROIZE  FOR  ACCWMING  S JM  00019100???? 

LDL(3)  ,KI  *  CRFATF  MOTFA  MODE  PATTERN  00019400???? 

LOA  prN.  %  TnRN  ON  PE-S  00019500???? 

IAL  *r3»  %  WITH  PEN  <  K  00019*00???? 

Sk-TC(3>  II  %  MnVF  I  BITS  TO  8C3  00019700???? 

CAND(3)  ,sVMnDE«i  AN!)  rtlTh  CURENT  MODE  PflTTFPN  00019800???? 

SfL(3)  ,^QnEA>  %  SAVE  In  MOpEA  00019900???? 

I  00020000???? 

LOOPJ*   LftLOl  .oFmoDE»*  PF-S  0  TO  M  ON  00020100???? 

LUEE1  S*3»  00020?00???? 

LUL(3)  .AD^AS  !S  ApnR  OF  MT<  (ppw  0)  OC02030O???? 

•CA0D(3)  *r2l  ■  *  AnOR  OF  J-TH  ROW  00020400???? 

UA  0(3)  I  *  J-TH  PO*  TO  RCA  00020500???? 

CSHL(3>  6|  11  FOPH  CU  AOOR  00020600???? 

LOL(O)  S*l|  %  I  00020700???? 

CSU8(0)  Sr;2i  *  I  -  J  00020*00???? 

CAD0C31  SrOI  *  AnOR  FF  (T-J)TH  ELFMFNT  Or  ROW  J  00020900???? 

LUA0(3)  SC3I  S  MnVE  fLEMrNT  L(J>I-J)  t0  SC3  00021000???? 

MLRN  $c3|  X  MULTIPLY  J-TH  Rf)W  BY  L  •*  J » I-JJ  0002U00???? 

I  00021200???? 

Li!L(3)  .sIXT^l  *  6a  00021300???? 

CSUB(3)  SCOj  %  6u    "    (I-J)  00021400???? 

RTl  $A»0(3)I55  RO'lTE  RESULT  LEFT  I-J  00021500???? 

lol(3)  .moofai  *  Turn  on  leftmost  off  pf  00021*00???? 

cshr(3>  1|  *  smipt  modfa  right  by  ove  00021700???? 

CUR(3)  .^OOFAJ  *  MfiVE  NEW  PATTFRN  00021800???? 

SI  L(  3  >  .^ODFAj  %  INTO  ►'ODEA  00021900???? 

LOEEl  Sr3j  %  AND  SET  E-RlTs  ACCORDINGLY  00022000???? 

I  00022100???? 

ida  s^l  *  rfaqy  for  sum  -  retrieve  00022?00???? 

adrn  tpmp|  *  roljtep  result  and  accum  00022300???? 

sta  tempj  *  sum  In  temp  00022400???? 

txltm(?)  »l00pj?  *  jump  if  j»sc2<i-1  00022500???? 

*  0O022600???? 
LDL(O)  ,SAVF0I  *  RESET  ACARO  00022700???? 
LDA  0(0)1  *  SUBTRACT  SUM  (\n    TEMP)  00022800???? 
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NEHROK 


SkSRN 
ST  A 
:LJL(2) 
LUEE1 
CALL 
LOS 

LlT(3) 
LOA 
DVRN 

su 

ST  A 

SI  IT(?) 
L  l'  A  D  (  2  1 
L  U  I.  C  2) 
C  A  k-  D  C  2  ) 
LoEEl 
LOA 
HlRN 
ST  A 

CIMPC( ?) 
LUEE1 
CLRAj 
SU 

jUMP 


tf^pi 

0(0)1 

t  pEO  j 

*C2f 

SyRTft'K 

Sa  I 

■1.01 
Sr3, 

*S  J 
0(0)1 

=  T  E  (•  &  I 

SC3? 

,  n  0  T  e  E  0 

.  S  V  ►«  r  u  E 
j>r2  » 
0(0)1 
*C3J 

0(0)1 


*  FROM  I-TH  ROW  AND 

*  STORE  bACK  IN  i-TM  ROW 

*  PATTERN  FOR  PEO 
X  GMI Y  P£0  ON 

U%    FINO  S&RT  OF  FIRST  ELEMENT 


S,*2 


%  TO  SAvE  TImE  LATER  cPE«TE  PFcTpROCAL 

t  OF  ROOT  (CAN  MULTIPLY  RATHFR  THAN 

%  n  t  v  t  n  f  w  H  r  n  f  n  p  m  i  n  r,  l  ) 

*  STORE  kECiPROCaL  OF  RO-T  jh    ROW  I 

*  Mp.VF  RESULT  OF  SORT  TO  SC3rl/LtI»lJ 

*  AmO  wjtH  HIRRFnT  MP'")E  3ATtFRn 

*  P  F - S  1  TO  M I N ( M •  N-I) 

<  LpAP  r<oA  ^TTM  (OLO  ROW  I  -  <;UM) 

*  Mi|L  T  hi     J/LtI»T] 

'(,  RFPLACE  A  RY  mfw  ROW  Or  L 

*  SAVE!)  MODE  PATTERN 
t  PF-S  M  +  l  TO  6  3 

3!  ZE^'O  FILL  REST  OF  ROW 


0(0)! 

MyTROW;  %      jlJMF  IF  I«*C1<N-1 


% 

% 

% 
% 
% 
% 

SOL 


A  IS  Nfjw  n'CTMansFD  TnTO  t  A r.  U  LT(L  TRANSPpSE'* 

THUS  >-E  HAVE   l  *  LT  =  F,  LET  ClT)X  =  y  ft  HO  THEN 

1,  SOLVE  LYSF  rOR  y  iFufc*ARO  SUHSTTTUT  TON) 

2.  SOlvE  (LT)x.Y  t-TR  X  (nACK     SU9STITIJTT0N) 


VEIT 


*FILL» 

StTE    E.OR 
CLRA: 
LDL(?) 
L0L(1 ) 
STLO  ) 
LDEE1 
LDA 

LOL(O) 
L  0  A  D  (  1  ) 
MLRN 
LDC(l) 

STOPrto) 

HTL 

LDL(l) 
CADO( 1 ) 


E  J    SETE1 .  F.  ANin.E  : 


.AnpAj 

.  PFOJ 

.SVmODE;* 

SCI  t 

0(2)* 

.AnRF» 

**C1  t 
*Cl  J 
*A| 

«.C1  I 

sa.i  ; 


xAnr:R   of    rEcOMpnsEp   mT* 

'A  MPpE  PATTFRN  FOR  PEO 
SAVE  TO  TnDICaTE  niRR 
TIRN  UN  oro  ONLY 

■1/L[0»01 

nCOsAnnR  DF  F  ARRAY 
%  C  1  s  F  C  0  ] 

Y[()]«FCO]/LCo.O] 
mOvE  TO  Af,AR  I 
pVcRVRlTr  F [ 0  J  WITH 


L 

:nt 


VOOE 


<*X*F> 


ROUTE  TO  pE  1  -  STORE  IN  TFMa  L^TER 


I 

FORLOUPt 


,svlod3j?;  loop  counter 
.omfi   %   loop  .frq^  1  to  n-1 


(S"1=PPW  INOEY) 


CADDC)   .nK'EJ    *  PUMP  ACAR  0  j  n 
CRTRTAC1)  cn«#FORMULT|  *  IF  I  GTR 


GET  NF^t  ELFMFNT  pF  F 
M  JUMP 


LDL(?> 

CSHRM) 

C0R(3> 


. SVMOOP; * 

If       * 

.SvnOoEix 


I  LEO  M 

^FT  CURRFNT  Mnr)E 
AK'D  ADO  I-Th  pE 


PATTERN 


00022 
00023 

0  0  023 

00023 
0002  3 
00023 
00023 
0  0  0  2  3 

00023 
00023 
00023 

0002a 

(10  0  2^ 

0002a 
on02a 
0002a 
0002a 
0002a 
0002a 

0002a 
0002a 
0002s 

00025 
00025 
0  0  0  ?  5 
00025 
00025 
00025 
00025 

0OC25 
0O025 
00026 
0  0  0  2  A 
000?6 
00026 
0  0  0  26 
00026 

nn0?6 
00026 
00026 

0  0  0  26 
00027 

OO027 
0  0  02  7 
C0027 
00027 
C  0  0  2  7 
POC27 
OP027 
00027 
00027 
0002A 
0  n  0  2  B 
0002B 
no02R 
0002B 
00028 


900???? 
000???? 
1 00???? 

?00???? 
300???? 

aoo???? 
50c???? 

60" ??7? 

700???? 
ftOO???? 
QOO???? 
000???? 

1 00???? 

200???? 
300???? 

aoo???? 

500???? 
600???? 
700???? 
ROC???? 
900???? 

000???? 
1 0O777? 

20077?? 
300???? 
/JOO???? 

500??7? 

600???? 

700???? 
POO???? 
900???? 
000???? 
10  0???? 

?oo???? 

300???? 

aoo???? 

500???? 
600???? 
700???? 
ftOO???? 

900???? 
000??7? 

1 00???? 
20^???? 
300???? 

aoo???? 
50c???? 

6  0  0???? 
700???? 
ftOO???? 
Q0  1???? 

000???? 
1 00???? 

200???? 
30'*???? 

aoo???? 

500???? 


3k 


LDEEl  *C3 I     *  PF-S  0  THRU  I  ON  00028600???? 

STK31  .SVMOO^It  S*VL  PATTERN  00028700???? 

iJR  TEMPI     3!  WITH  PE  I  ON  jtR  CAN  N"W  RE  STOREn  00028800???? 

%  TF.MPB(0»V[I-n»Yri-P3»..Y[ol»0...0>  00028900???? 

CLC(2>J  0O029000???? 

CaDP(?>  sell  00029100???? 

LDA  *C?I      %  ROA«I  00029200???? 

SBM  peni     *  RGA«|-PEM=(IiI-l»..«01  00029300???? 

SkN  AP  t  OP0294OO???? 

LDX  <B(      %  MOVE  TO  RGB  FOR  RGY  00029^00???? 

*  000296OP???? 

EORMULTj  LDL(3)  ,AOPA|  %  tC3sAi)DR  OF  L  MTX  00029700???? 

i-DA  *0(3)l    *  °Ic!<  OP  "OW  I  OF  I  00029800???? 

LDS  kAj      %  SAVE  LlI.Tl  IN  $S  FOP  IATFR  USE  00029900???? 

*LRN  TEMPJ    %  COKPUTF  LrI*Kl*YrK]   *»I-f »...I-1  00030000???? 

I  00030100???? 

C  L  C  C  3  ) ;  SLIT(3^  sROKSMl  EXCHL(3)  SICRl  00030200???? 

I  00030300???? 

LDLCZM  .PEOl    %  MOPE  PATTERN  FOR  PEO  00030400???? 

LDEE1  *C?r      *  TURN  ON  PFO  00030500???? 

LOAD(O)  <C?I      ?  *C2sFCI]  00030600???? 

5>BR^  ^C?  *      »  !>ri&    OF  PTO  a  FCI3  -  l<"I»Kl*FtK)  00030700???? 

CHSA;  *  WHFRF  K«I-M. ...  1-1  00030*00???? 

^LRN  sSl       *  KlJLT  AHOVE  HY  1/Lrl.P  SAVFD  00030900???? 

'/  I-N  SS  PRTOP  TO  LHG  SUM  00031000???? 

LnC(?l  *A|       %  MOVE  in  ACAR  2  00031100???? 

STORFfU)  eC?l     %  OVERBITE  Fill  WITH  YM]  00031200???? 

LDLOl  .PFO?  00031300???? 

CnMPCf3)j  *  PF-S  1  TO  63  00031400???? 

inEEi  %r^:    %  turn  on  all  pe-s  EvcE3T  peo  ooo3i5oo???? 

LOA  TE^PJ    *  PRINT-  TEMP  TO  RGA'  00031*00???? 

L0L(3>  .SVMUDrJ*  TURN  QN  5AVE0  00031700???? 

,  LDEE1  uC3i     *  MPOF.  PATTERN  00031800???? 

kTL  SA.1I     7  "r'ft«{  »I  I  J»  Y[I-1]»  ...Vri-M7,0)  00031900???? 

LESSTA(l)  SpH.B'IMPlj  %     IF  I<1  WAIT  TO  ST^RE  F,r,R  00032000???? 

STR  TEMP.I    *  I  GEO  K,  TEMP  CONTAIN5  LAST  M  YtK]  00032100???? 

IDA  %Xt                %  RU^P  RGX  BY  ONE  00032200???? 

LDL(3>  .ONE?    *  TO  Pi CK  UP  NE*T  ROW  00032300???? 

AqH  5C3t     ?  nr  L  00032400???? 

SWAP;  %  GO  THRU  RGB  AS  NO  00032500???? 

LDX  kBj      ?  DIRECT  PATH  FOR  LDX  $fi  00032600???? 

BUMPU    TXLTA^(l)  .Fn^L^OPJ  00032700???? 

I  0003280O???? 

%  BACK  SUBSTITUTION  "  SnLVE  LT  *  X  c  Y  FOR  X  00032900???? 

I  0003300O???? 

St-TE  E.O*.-El  SEtE1  E. AND. El  00033100???? 

CLRAj  00033200???? 

LDLC2)  .AnRAj   *AprR  CF  pFCOMPOSEo  MTX  LT  00033300???? 

CADD(?)  .N1>      *  LAST  RDw  OF  LT  00033400???? 

LDL(l)  ,PFO»    %  MODE-PATTERN  FOR  PEO  00033500???? 

STL(l)  .SVMODEjX  SAVE  TO  INDICATE  CURRENT  HOOE  00033600???? 

LDEEl  %<;\t             *  TURN  UN  pEO  ONLY  C0033700???? 

LDA  D(2)J    *  =l/LtN"l»N-l ]  00033800???? 

LDL(O)  .ArPFl   %  JCOsADPP  OF  Y  ARRAY  (F  DESTROYED)  0O033900???? 

CADDfO)  .Nil      J  AOOR  OF  Y[N-1J  00034000???? 

LOAO(T)  Ull      *  *f.l  =  YlN-1]  0O034100???? 

^LHN  *Cll     %  Xf N-l ]"YfN-l 1/LtN-l .N-t 1  0O03420O???? 
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Lncci > 

STORE ( 0) 
KTL 

LDLd  ) 

cAunf i ) 


SAj  t    MOVE     TO    ACAR    1 

tCl>  *    OVERKRlTf    YtN-11     WITH    XCN-ll 

XA,U  S    POtiTt    TO    PE    1    -    STORE    IN    TFMP 


LATER 


.SVLC00;*    LOOP    COUNTER 
,ONF?  *    L^OP    COUNTER    e    1     Tn    N-l    Bl 

%     ACTUALLY    GO    FROM    N-?     T0    0 


T    WILL 


BKLOOP{ 


CSUPfO)        .OK-Fl  *    OECR    ACAR    0    TO    GET    NFVT    FLTMFNT    OF    Y 

fcRTPTA(l)     $n«>B"MMLT|        *     IF    *C1>M    Jliup     (T,E,     T<N-M-1) 

%     SINrF  I-(N-1 )-$*l,   *C1=N"T-1 
*  So  *Cl>M  MFANS  "•I-1>M  Or  I<N"m-1 


BKHUl.Tt 


LPL(3> 

CSHKO) 

C  0  R  (  3  ^ 

LDEE1 

STL(?l 

STR 


L0L<3) 

CADD(  3) 
CSUP(3) 

LnA 

LDS 
>'LRN 


,svMODr:* 

i  t  % 

,SVMODpl* 

tEn-p;   7 


I  LEO  M 

GET  CURRENT  "nDE  PATTERN 

ANqAjdT-ThpE 

PE-S  0  TMRU  I  ON 

S  A  V F  PATTERN 

MTH  PE  I  ON  $R  CAN  N^-/  RF  STOREo 

TF«Psco,y[N-T  j»xrN-T  +  n,.yrN-i  3...0) 


.  appA  ; 

.M  : 
tci ; 
o(3W 
*Aj 

temp; 


fC3«AUl)R  qF  Lt  'U" 
PICK  UP  ROW  I      I=N 
COPRSP  VALUE  OP  SCls 
pn  A  =  Riiw  I  OF  LT 


1»   ?, 


1.   0 
N-2»N-1 


%    SAVE  Ltl.I]  IN  $S  FOR  LATPp  USE 

%    COMPUTE  LTCI»K]*y[Kl   Ksl  +  1,,,,UM 


CAl 


C H LATE    RnwsUM 
CLC(3ii       SLTT131    sROWSMJ    MCHLO)    SICRJ 


LDL(?)  .PFO;  K 

LDEF1  »C?»  * 

LO  A P  f  0  J       j C  ?  :  *  ■ 

C  h  S  A  J  z 

►LRN  %S!  * 

LDC(?i         FA;  * 

STDhE(O)    tc?f  * 

ldlo>        .Pro; 

tOMPC(3)j  * 

LUEFI  *r3,  5r 

LDA  TEf'PJ  * 

lpl<3>       .svmod"!* 

L0E.E1  %C3t  f 

K  T  L  j  a  ,  1 J  T 

LFSSTA (J  )     fnfl'*1  J 
S  T  R  T  E  >'  P  I  * 

TXLTA'iCl)     .PkLO^PJ 


MOOE    PATTERN    FnR    PEO 

TORN    JN    PEO 

*C?=Y[ I ] 

pr,A     D  F     pFO     =     Y[Il     -LTri,K1*XrK] 

W  H  E  R  F     K  =  T  +  1»...I«-M 

Ml'LT    ABOVE    pY    l/l  r  I  »n    SAVFD 

IN    iS    PRIOR    TO    LfJGSUM 

MOVE    TO    ACAR     2 

OVERWRITE  Y  T I  ]  WITH  X  r I ] 


PF-5  1  TO  63 

TURN  ON  fttL  PF-S  FypE°T  pFO 

WRING  TE^PTO  RGA 

ThRN  DN  SAVEO 

MPDE  PATTERN 

pgas(Xct  i,xci+n»      xrr+M3#o) 

%    If  I<M  WAIT  TO  STTpE  pr,R 
I  GEG  m.  TEMP  CONTAIN*  LAST  M 


xt«i 


LDLC3)  ,,|SR"DE;^  rfstore  usfr  mooe 

LUEE!  $-3» 

LL'LC(M  . uSp^O;  %    RFSET  ACAR  0 

LlLC  1  i  ,uSr<*1  I     *  RESET  ATAR  1 

LLL(3i  .RETliBj  *  PICK  iiP  RFTURN  A009ESS 

EJCMK3)  Ficpj    %    AND  RETURN  TO  OUTSIOE 

FN0  h.sdi  vEfti. 


P003A 
0003A 
0003^ 
00034 
00034 
00034 
00034 
0n03«> 
0003S 
0003S 

0^035 
00035 
00035 
00035 
00035 
00035 
00035 
00036 
00036 
00036 
00036 
00036 
00036 
00036 
00036 
00036 
00036 
00037 
00037 

0  0  0  3  7 
00037 
00037 
00037 

00037 
00037 
O0037 
00037 
00038 
00038 
00038 

00038 
00038 
00038 
0  0  0  38 
00038 
00038 
00038 
00039 
00039 
0003? 
00039 
00039 
00039 
00039 
00039 
00039 
00039 


300???? 
400???? 
500???? 
600??7? 
700???? 
POO???? 
900???? 
0  0  0  ?  ?  ?  ? 
1 00???? 
200???? 

3  0  0  ?  7  ?  7 
400???? 

5  n  o  7  7  7  ? 
600???? 
700???? 
8 00 ? ?7? 

900???? 
0  0  0???? 
100???? 
?00???? 
30<v???? 
400???? 
500???? 

6  0  0???? 
700??77 
*00???? 
900???? 
00^777? 
1 O0??7? 
700???? 
300???? 

4  0  o  7  7  7  ? 
500??7? 
600???? 
700???? 
SO!.)???? 
900???? 
000???? 

1007??? 

200???? 
300???? 
4  0  0???? 
500???? 
6  0  0  7  7  7  7 

700???? 
600???? 

90n??7? 
OOo???? 

100???? 

200???? 

300???? 
400???? 
500???? 
600???? 
700???? 

flO',???? 

coo???? 


00040000???? 
OOOftOlOO???? 
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APPENDIX  B  -  BS0LVEL28  ASK  LISTING 


?U$ER»CACIQUS  OOOOOlOO???? 

7C0HPILE  ►lC0/A^K3/rtSnLvE128     KITH  ASK  LIMAKY  00O002O0???? 

?BCL                 00000300???? 

fct-GlN  oooooaoo???? 

CHhS      641             '''  00000500???? 

FILL      l?8i              •  00000600???? 

*  DEFINE  PRINTAb  00000700???? 
%  SETH  F.Or.E,  SFTj  El.ns.El;  SETE  E.OR.-F*  SETE1  E.ANn.Fi  O0000fl00???? 
X  STA  TFMPT;  STL(3>  .sAVr3j   LjT(3)  *1 • TfMPT+63 *  T FMPT I  DtSPLaYR  *C3»16;  00000900???? 

*  CU<3)>  LDEE1  *C3*  I ULf 3  )  .SAVE3J  SETE  H.OR.E?  SETE1  J.OR,El*#t»  00001000???? 
ZEPHl  EUU  tnOi  t  ZrRQ  ooooiioo???? 
ONEi  EUu  $r>lj  *  FTxEn  POINT  ONE  00001200???? 
SIXT4:  EUU  St)2i  X  FIXED  POINT  64  00001300???? 
TWO»  EyU  Sr)3l  *  FIXED  POINT  2  00001400???? 
ADPAt  L«U  Sn4|  *  ADDRESS  OF  HTX  00001500???? 
ApPF*  Et»U  S05|  *  AnriRtSS  of  rhs  oooouoo???? 
LASTMJ  EUu       *q6>     X  INDEX  (jF  f.\|-M>TH  ROW  OF  MTY  00001700???? 

*  I.E.  FIRST  OF  LAST  H  ROWS  P0001800???? 
K|      E^U       Sr>7,     *  USED  IN  C&FATING  MODE  DATTFRNS  00001900???? 

*  NOTE;  1  LEO  K  LEO  M  00002000???? 
Ml  EUU  $n8t  %  BANP  *IfHu  (  INPUT  )  00002100???? 
N|  E«U  Sr>9|  *  DIMENSION  OF  MATRIX  00002200???? 
Nlj  EUu  S',101  *  N-l  00002300???? 
hOliEAl  EUU  *i)13j  *  40PE  FATTFRN  FflR  ROw  1  00002400???? 
MODEBJ  EUu  SnU?  *  "nOE  PATTERN  FOR  ROw  1*2  00002500???? 
PEO|  EUIJ  $)15J  X  MnOF  PATTfPN  FfiR  PFO  OvLY  .  00002600???? 
MOf.Ell  Et»u       *^16»    %    MOnE  pATTFRN  FDR  FIRST  64  pR  ,  ESS              00002700???? 

%    ELEMFMTS  np  a  GIVEN  RP>'  00002*00???? 

M0DE21  EUu       S o  1 7  J    *  MODE  PATTFRN  FOR  LAST  v,-6*i  ELFMENj-S  00002900???? 

*  Or  a  GIVtm  row  000030O0???? 
NOTPEOtEUu  S-)1R»  %  r-'ODF  PATTERN  rpR  PF-S  1  TO  ft3  00O031OO???? 
SVMODE»EUU  S'UPj  *  SAVE  MODE  00003200???? 
SAvEOj  EUu  S^20;  *  SAVE  AcARO  C000330O7??7 
SAVE3»  EUU  Sn22t  *  SAVE  ACAR3  00003400???? 
SA^ADPtEUu  $T23:  *  AO^  TE^P  00003500???? 
SVLDOP{EWU  S12&I  <  SA^E  LOOP  INDEX  (0  TO  *'-l)  00003600???? 
RETuRj  EUu  S->25;  X  RETyRN  AqoRESS  00C0370O???? 
USRMDEjEUu  Si>2f.,  *  IJSFR  KQDE  PATTERN  00003800???? 
SiXT3l  EUu  S>)27«  *  FTX-PT  63  00003900???? 
ADRAN1,EUU  t)2Bl  %  INTX  OF  LAST  RpW  OF  A/L  M.TX  00004000???? 
USRC°»  EUU  S.i29|  X  'iSER  AqAR  0  00004100???? 
USRClJ  EUU       *r>30J    X  USER  AGAR  1  00004200???? 

*  00004300???? 
TEMPI  PLK  li  %  TEMPORARY  ROW  STORAGE  00004400???? 
TEMPI  »  BLK  1j  %  SECOND  ROW  OF  TEMP  STORAGE  00004500???? 
TEMPRJ  BLK  lj  X  AUXILIARY  REGISTER  P0004frO0???? 
TEMPPJ   BLK       1,        *  TF^P  STORAGE  00004700???? 

*  TE^PTiBLK  1,  X  FqR  nE8  IGGING  ONL*  00004800???? 
RGXU  BLK  l,  |  TEMP  REG  FoP  RGX  00004900???? 
RGX?|  PLK  1,  %  TpMP  PEr'  ?  FOR  RGX  00005000???? 
PEN:     CATA  0. l»2.3»4»5»6»7,fl.9. lo. J 1# 12. 13# 14. 1?» 16# 17.18»  00005100???? 

5('#?T,?1.22»23»?4,25»26.27.?fl,29.3*.,l»32.33#34,  00  005200???? 

?S»36.37,38»39.4r,  ,A)»«?.43.44.tt5»46.f'7,4fl,49»50,  on  0  05300???? 

5<»52.53»54.5S»56»57.58.59»6r>»6l»62>*3J  000  05400???? 

ft  0000550C???? 

*  SUBROUTINF  rUwSm  caLC'LATES  T»-E  Shm  Of  THE  ELEMENTS  00005600???? 
ft   IN  RGA.  RGA  MLl  CONTAIN  THE  RESULT  OF  THE  SUM  UPON  00005700???? 
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%    EXITING    THE     RO.iTlNF 
t 

KOMSMt  > 


.SAVE3» 


FILLI 

STL(3) 

CUMPC(?)i 

LDEEl 

TLRA  j 

SLTE 

SLTEl 

SLIT(3> 

HT| 

At.RN 

CSMLC3) 

T  *  L  T  A  M  (  2  >     ,  -  4  t 


E  .  A  N  D .  r  t 

1»6»  1 J 

1  ? 

*A,r>(  31  » 


%  SAVE'    RETURN    AOOR 

?.  PF.VFwSE    MODE    PATTERN    IN    fc2 

%  Tt.iRN    ON    pE-S    NOT    nESTREr>    rflR    SUM 

%  AM)  CLEAR  THTIR  RGA 

x  turn  on  all  pe-s 

s  for  logsum 

%  LOOP  COUNT 

*  SHirT  imoex 

*  ROUTE  AMO 

■?  a  I',  n 

?  DnUPLE  SHIFT  OIST 

%  CONTINUE 


L0LC3) 
EXCHL( 3) 


,SAVE3» 

*Icr; 


GET 

AND 


RETURN 
KETURN 


AOOR 


I 

*   ENU  pE  SUbRiUTI^E  POWSM 

I 

BSOLVElpfclENTRY) » *       FTLL« 

S'L(3)  ,»ET"P»     *    SAVE     RETURN    A">PRES$ 


S  E  T  C  (  3  )       El 


*    GrT    USfP    MODE    PATTrRN 


SU(3> 
S  T  L  C  0  ) 
SlL(l) 
SETE 


•  uSP^'DE  '  * 
.iiSprOj  % 
.  ijSPrl  ;     % 

E.no.-f? 


A  »■•  n  savE  TT 
SAVE  ACAR  0 
SAVE  AcAR  t 
S  f  T  E  1   F. .  A  N  n  .  r  ; 


CLC(0>: 

STKCh   .ZERn> 

cumpc( m i 


S  1  L  (  0  ) 

LlT(O) 

ST  L  <  0  1 

LlT(C) 

SlL(O) 

LlT(O) 

SlL(O) 

CSUB(0> 

ST  i.(O) 

LiT(O) 

STLCOl 


,OnFl;  *  INITIALLY  ALL  PE-S  ON 


■1  5 

.ONE' 

a?t 

•  TaO| 

=  <4  J 

.SIXT4* 

•  n  n  e  i 

.Sl*T3| 

s  1  u  0  0  0  0 
.PEO  I 


%    FlyEo-pT  TwO 


%    FIXEO-PT  63 

ooooooooooooooomi 
*  i10de  for  1  pe 


00005 
00005 
00006 
C0006 
00006 
00006 
00006 
00006 
00006 
00006 
00006 

00006 
00007 
00007 
0  0  007 
00007 
O0007 
00007 
00007 
O0007 
00007 
OP007 
00  0  0  8 
00008 
00008 
00008 
00008 
0  0  0  08 

noooa 

00008 
00008 
00008 
00009 
0  0  00  9 
0O009 
00009 
0  0  0  0  9 
0  0  00  9 

00009 
00009 
00009 
00009 
00010 
00010 
00010 

onoi  o 

00010 

oooio 

00010 
00010 
00010 

oooio 

00011 
00011 
00011 
00011 
00011 


*po???? 

900???? 
0  0  0???? 
100???? 
200???? 

3  0  0???? 

400???? 
500???? 

6  0  0???? 

7on???? 

800???? 
900???? 

0  0  0???? 
10*???? 

200???? 

300???? 
a  o ',???? 

500???? 
600???? 
700???? 

800???? 
900???? 

ooo???? 

1 0  0???? 
200???? 
300???? 
400???? 
SO'1???? 
600???? 
700???? 
800???? 
900???? 

ooo???? 

1  o<>???? 

?oo???? 

300???? 
400???? 
5  0  0  7  ?  ?  ? 

*00???7 
700???? 
8  00???? 
90(i???? 
000???? 
IOC???? 
200???? 
300???? 
400???? 
500???? 
600???? 
700???? 
800???? 
900???? 
0  0  0  ?  ?  7  ? 
100???? 
200???? 
300???? 
400???? 
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CSHR(O) 

ST|.(0> 

ALIT(2) 

l.UA(j(?) 

STL(O) 
ALIK2) 
LUAlK?  } 
LUAO(O) 
STL(O) 

CSljO(O) 

STL(O) 

ALIT(2> 

Ll»AD(2) 

LUAO(Ol 

STL(O) 

LDL(O) 
CbU3(o) 

SIL(O) 


$COl 
6| 

•  aora * 

ll 

$co» 

.ADPr; 
It 

$C0  I 
SCO; 
.Mt 

•  nNFI 
.  Ml  i 
U 
*CX)| 

.^t 

•  ,; 
.1; 

.lastm 


LUL(O)  .mj  it 

csus(o)  .sIxt4;  2 

L^A  PpNj  * 

I A G  SCO :  * 

StTc(O)  It  ? 

CUMPCO)  I  « 

LUEEl  ScO:  * 

S^LCoi  »m[iof2;  x 

STL(O)  .sVhnoEf? 


LUL(O)    .oFot 

CUMPco): 

STL(O) 

LlTCl  ) 
CADDd  ^ 
CHOTUl) 
STLCll 


•NOTpEOj* 
0,1.05  % 
•  M  It      * 

•SVi  nCPi 


ALIT(?) 

LUAD(2> 
LUADd  ) 


It 

srl  t 


M  ADR 

NEXT  PA'AMFTCR  aDDr 


Annp  pf  MTX 

SHIFT  CU  ADDR  TO  GFT  P^W 

pf  AnnR  nF  mtx 

Hi.iMP    LIST    POINTER 

AnDR    OF    F    VEcTOR 

STORE    CD    «nnp    jnj    ADR 

RllNP    ACAR?    FOP 

PIC*  HP  PARAMETER  in  $ro 

value  in  AnoR 

DIMENSION  OF  MTX 


H-.l 

Rl)MP  ACAP^  FOP  NEXT  PAOAMFTFR 

PICK  LP  PARAMETER  IN  tr0 

VALUE  IN  AOOR 

BAND  WiOTH 


AODR 


N-M 
IMDEX 


UF  fM-M)TH  ROW  Or  Mtv 


3A>'D  wIDTH 
M-64 

°E  INDlCirS 

1   fits  on  *hfre  pfm>ravid  wjptm 

MOVF  I  BHS  TO  ACARO 
RFVERSE  PATTERN 
PF-S  0  TH  m  -  64  ARF  0« 
Sa'vL  in  MpOE2 

also  save  in  svmodf 


PF-S  1-63  ON 

SAVE  IN  NlTpEO 

SrT  UP  LOOP  COUNTER 

To  GO  FRO*  0 

IJP  TO  MM  (BY  ONE) 


%    Annp  [F  FLAG 
%    CONTENTS  NOW 


PARAMETER 
IN  ACAR  1 


ZEpTA(l)  »+l)     %    IF  FLAG=0  THEn  DECnMPncE  MTX 
JUMP      SOLVFIP*  OTHERWISE  JUMP 


BEGIN  MAIM  LOOo  .  REOUCE  A  TP  LO"EP  TRIANGULA* 

-  *C0  CONTAINS  AOORESS  OF  ROW  TO 
BE  R£pLACEL  (A  IS  DeSTPOVED^ 

-  tr1  IS  Lf'Up  COUNTER  .0.1, ...N-l 

SETE  E.OR.-Fj  SE^El  F.ANp.El 

LlT(l)    0,1, Ot   <    SET  UP  ACAR  1  -  COPRSP.  TO 

CAD,(1,   .ml  |     %    TO  L.OCP  FPOM  0 

CHOTL(l)  2d |      X  UP  TO  N-l  (BY  ONE) 


C0011 

noon 

00011 
oooi  1 
00011 

00012 
00012 
00012 
00012 
00012 
00012 
00012 
00012 
00012 
00012 
00013 
00013 
00013 
00013 
0O013 
00013 
00013 
00013 

00013 
00013 
00014 
00014 

00014 
00014 
00014 
00014 

0001* 
00014 
00014 
00014 

00015 
00015 
00015 
00015 

00015 
00015 
00015 

00015 
00015 
00015 
00016 
00016 
00016 
00016 
00016 
00016 
00016 
OOOlA 
00016 
00016 
00017 
00017 


500???? 
600???? 

700???? 

Ron???? 

900???? 
0  00???? 
100???? 
200???? 
300???? 
400???? 
500???? 
600???? 
700???? 
800???? 
90o???? 
000???? 
10U???? 
200???? 
300???? 
400???? 
500???? 
ft 00???? 
700???? 
RO'J???? 
900???? 
000???? 
100???? 

200???? 
300???? 
400???? 

5  0  0???? 
iSOO???? 
700???? 
POO???? 
900???? 

ooo???? 

100???? 
200???? 
300???? 
400???? 
500???? 
600???? 
700???? 
flOO???? 
900???? 
000???? 
100???? 

200???? 

300???? 

aou???? 

500???? 

6  00???? 

700???? 
flOO???? 
900???? 
000???? 
100???? 


39 


L  L>L  t  1  5  .SVLHOPIX  SAVE    FOP    USE    !M    FOR/BAf*K    Sl'pST    LOGIC 
I 

lul(O)  .appaj  *  first  rcn  is  spec  case 

IL>A  0(0)!    *  LOAD  IT  IW  PfiA 

JIMP  NfWpnKi  X  JhPP  PvfR  MULTIPLYING  LOGIf.,  ETc» 
* 

NXTRn*5!LLL(0)  .AOPA?   *  AnJhST  ACARO  TO  CPNTAI«  TwF 

rAOO(n)  $ri|    x  ahtress  "nr  the  i-tm  Rn*1  of  mtv 

0*00(0)  J,rl|             %  2    FL'rtS    Uf    PE"M    pER    RPW     nE     A/I 

LDL(2)  .nNE!         *  SET    k  =  1.       wRL    ONLY    fh«NGE 

SU(2)  .k;             *  "HEN    l<6<» 

U'L(3)  Scl  S            %  I 

ffcUHO)  .mj      X  I-f-  IN  fc3 

LESSTA(l)  io^r^Ml  ?  JUKP  IE  I   <   (N-M) 
I 
I   EOH  ThE  "LAST  M  pOkS  C'EATF 

x  specific  Mpnt  patterns 
i 

lllc?)  .•,nnr2;  *  i  r.EQ  n-m 

2tPTA(?)  .s-TrTl:     %    jump  Tf  MPrE2  FALSE  (T  r,£Q  N-M) 

r.ShL(2)  li       %  Tur.N  rF P  PE    N-i-6« 

ST  L ( ?  ^  ,-opr2;    I  Awn    SAVE    PATTEpN    IN    mp*F2 

S^lP  .C«Af,3j 
I 

SHIFT1*    L.DLC2)  .,-nnEl?     X  I    rF0    *-6a 

C.ShL(?^  Ij                  i  TnPN    PFF     PE    N-I 

SU(?i  .-pnrlj    *  Ami:    S/wE    pattern    in   Mpr>E    i 
X 

CKAC3«       LtSSFM)  S^6,rIXAC?JX    J'JHP    IF    M    LFQ    I     *    I    <    fM-M) 
X 

X  FpP    FIRST    M    9f)Vs    '    A "»  itJST    VALUE     OF    K    AND    SET 

X  L^'FR    LIMIT     pr    LPPP    aO    (PP'1    0) 
X 

LLSSF(t)  i^2,clZ*R3'>     *    JUMP    T  F    fi/J    LEO    I    <    M    {K-n 

L C: L  (  3  ->  .  sI*t4>     %  64 

CS|;P(3)  frMj            *  6<)"I 

STlO)  ,k>               %  K    e64    -    I 

CLEARS*    CU131I  X  ST/-KT    AT    FIRST    R0> 
X 

FIXAC2J :LlT(2)  0,1, rJ       X  SET    UP    Af.AR    2    «    CPPPSP.    Tp    j 

C  A  D  D  (  ?  i  Sell            *  Tn    l.ppP    TP 

CS|.B(2)  .PNF»         X  1-5     (SCl-1)               (UPPEp    lImTT) 

CHPTLC?)  2a;               X  FRO*    MX(0,IrM)     (LPwEP    LlMjT) 

OPD<2>  Sr3J 


THE    FLLLPv-I^n    i_OPe    IS    TRAVERSED    m-  TI'«ES 
EXCFPT     rtHEM    I<M    -    TmiN    ONLY     I     TI«ES 

STL(0i  •  c; A \/ F 0 J     %    SAVE     ACaR    0 

StTE    F..0R.-EI    SETE1     E.ANp.El 

CLPAl 

STA  TEMPI  %    ZFFPI7E    FOR    ACCUMING    S'.iM 

STA     Tempi >  *  also  zepptze  secpnd  rp^ 

LCL(3)  .Kl 

LbA  PfM 

IAL  sc'3|     X  I  PITS  SET  WHERE  PFN  <  K 

StTCO)  I|       X  MHVE  TU  ACAR  3 


0001720 
0001730 
0001740 

oroi 750 

000176,0 

0001770 
00017*0 
0001790 
0001800 
0  0  0  1  P  1  0 
0P01fi?0 
0  0  0  1  fi  3  0 
0001840 
0001850 
0O01860 
0001870 

nr>oi8AO 

0001890 
0001900 
0  0  U  1  9  1  0 
0001920 
0  0  019  3  0 
00019£0 
0001950 
00019*0 
0001970 
00019RO 
0001990 
0  0  0  2  00  0 
0  0  0  2  010 
0  0  0  2  0  ?  0 

00  020  30 
0  0  0  2  0  4  0 
OO02O50 

no020f 0 
0002070 

00020 fO 
0  C ■  0  2  0  9  0 
A0021 00 
0O0211O 
0002120 
0002130 
0  n  0  2  1  4  0 
o  n  0  2  1  5  0 
0002160 
00  02170 
0002180 
0002190 
0002200 
000221 0 
0  0  0  2  22  0 
0  r  0  2  2  3  0 

0002240 
0002250 
C  r,  0  2  2  6  0 
0002270 
00022eO 


0???? 
0???? 

o???? 
u???? 

0???7 

()???? 

0???? 
0???? 

r???? 

0???? 

0???? 
0???? 

( '•???? 
'  ■•???? 

0???? 

(!???? 

()???? 

0  v  ?  ?  ? 
I)???? 
r  ?  ?  ?  ? 

l»???? 

u???? 
('???? 
o???? 

0???? 
n???? 
0???? 

"???? 
o???? 

0???? 
0???? 
0???? 
07??? 
"  ?  ?  ?  ? 
0???? 
0???? 

()???? 

0???? 

0  7  77? 

O7077 
0???? 

0???? 
0???? 

0  7777" 
07??? 
0???? 

0???? 

07777 
(■???? 
(???? 

0???? 

0???? 

0???? 
0???? 

0777? 


UO 


C*KD(?)   .uonrli  *  Awr  with  euRREK'T 

STL(3)    .MODEAl  %    SAVE.  IN  MQDEA 


LDPPJ**  StTE  F.&R.  -F;  SETE1  F. AND. El 


% 

LOoPl » 


LOLC3) 
c*pp(3) 

CAOUC3) 

LUL(O) 

CSUB(O) 

L()A 

Li>S 

CSHLC3) 

CACU(3> 

LL'ADOl 

LESSTO) 

C&OBCO) 

Lt'A 

M|_PN 

LK(3) 

Lf;A 

L  0  L  C  3 1 
CSHL(O) 
L «-  E  E 1 

A  f  R  N 

STA 

JUMP 

MLRN 
S1A 

LUA 

MLRN 

LOS 

LDLC3) 

CSU8(3> 

LUA 

RTL 

LULC3) 
LDFEl 

LUA 

AOPN 

STA 

CSHR(3) 

CUR(3l 

S1LC3) 

LULO) 
CRHTLO) 
C  A  m  D  (  3 1 
Si L ( 3) 

LUI 
StTE 


•  aOra 
sol  J 

*C2l 

f»(3), 
1(3)! 

*r0, 

SC3J 


Anr.R 
Anr.P 


rF 

PF 


mtx  (Rpw 

J-TH  ROW 


MODE  PATTERN 


0) 


Twf«  RP*S  OF  PEM  PEP  RO''  Of  a/l 

I 

I-J 

J-Th 
RFST 

fprn 

appr  pf  (t-j)th  element  of  row 

HpVE  ELEMFNT  L(J»I-J)  T0  Sf3 


POW  TO  PGA 

DF  ROW  TO  RGS 

fu  AOOR 

PF  ( T-J)TH  ELEMENT 

ELEMFNT  L(J»I-J)  T0 


*02,l  PPD1*  "*  Jl'KH  IF  I-J  <  64 


.sixt** 

Sr3: 

.  S I  X  T  4  I 

$  c  o  J 

SAtO(3)  »  » 


TfMp  j 

.qvrnDE 
0(0)  J 
$03; 
$RJ 
TEMP? 

EnDJ* 

Sc3 , 

Tfmpp: 

4s* 

tr.3» 

*a; 


ir.0*  * 
TfMPp,  % 
SA»T(  3)  t  g 

•  -40RFA.)  '* 
sr3t     * 


*R» 

T  E  K  P  , 
TrMPt 
ll 

•  MOnr'Ai 
.^nnFA; 


•SVunPEt  t 
0(0)*    * 

•Monrii  * 

•  MOrFB;  * 
V3,  X 
I .  ANn. -? i 


trfBf  r.U"64 

RFST  f'F  RHW  TO  RGA 

MiiLT  P.f  LU»I-J) 

Rru.TE  Kt'SULT 

I-J  LEFT  INTO  *R 

FROM  ir  A 

ACfUM  SUM  IN  TEMP 

SHIFT  jjRlR  M0PE2 

I-..'  LEFT  END  pFE 

CpFATF  ne-  mope 


1?  STrfc-E  bAC<  IN  TEMP 


niLT  J-TH  RPW  ry  LtJ»I"J] 
SAvF  IN  A|Jx  REr, 
SrCflNI''  PART  PF  ROW  J 
SECPNn  PART  PF  RPW  J 
SaVe  FnR  LATER  USE 


*   LU.I-J] 


64 

6«  -  CI-J) 

first  part  pe  row 

RnijTE    RESULT    PF    FIRST    "ART    PF    POW 

*    L(J»1"J) 

T|  PAD  &C3  WITH  MOOFA  PTRN 
AND  SET  E-BITS  ACCORDINGLY 

READY  FOR  SUM  -  RETRIEVE 
RPUTEn  RESULT'  AMD  Acmv 
SUM  I'M  7FMP 

TljRN  P\  LFTTMOST  PFF  PE 
PP  WITH  PLD  MPnEA  AND 
AMP  SAVE  TN  MpoEA 

MnPF  PATTERN  FnP  LAST  M"6A  FLMTS 

shift  i-j  left  enp  arp'nd 

AMP  WITH  niRRFNT  MppE  °ATTFRN 
SAVE  In  MPnEB 
LOAD  I  PITS  FROM  ACAR3 
SFTE1  E. AND. El 


0002? 
00023 

0  0023 
00023 

00023 
00023 

00023 
0O023 
00023 
00023 
00023 
00024 
00024 
00024 
00024 
0O024 
00024 
00024 
00024 
0  0  0  2  4 
00024 
00025 
00025 

On025 
00025 
00025 
00025 
00025 

00025 
00025 
00025 
00026 
00026 
00026 
00026 
00026 
00026 

00026 
00026 

00026 

C0C26 

00027 
00027 
00027 
00027 
00027 
00027 
00027 
00027 
00027 
00027 
00028 
00028 
00028 
00028 
00028 
00028 


900???? 
000???? 

JOO???? 
200???? 

300???? 
400???? 

500???? 
600???? 
700???? 
R00??7? 
POO???? 
000???? 
100???? 

?oo???? 

300???? 
400???? 

soo???? 

600???? 
700???? 

POC???? 
90C???? 
000???? 
1 00???? 
200???? 
300???? 
400???? 
500???? 
600???? 
700???? 
POO???? 
900???? 
000???? 

100???? 
200???? 
300???? 
400???? 
SOO???? 

600???? 
700???? 
800???? 

Qoo???? 

000???? 
100???? 

2007777 

300???? 
400???? 

500???? 
6O0???? 

700???? 
BOO???? 

900???? 

000???? 
100???? 

200???? 
3007777 

400???? 
500???? 


1+1 


ENDJ» * 


NEWROM 


CSu«( 3) 

rtl 
loa 

ADRN 
Si  A 

LDL(3) 
I  i  P  T  A  (  3  ) 
C*M)(31 
CAf.D(3) 
LLTE1 

A[,RN 
STA 


•SlXT^I     *  SHIFT     *A    LETT     I-J 

*c.O: 

$Si()(3)I31  RnitTEn    DATA     IN    $R 

T E M P  J  6  A n D    Tn    TEMP 

rr; 

Tempi        *  save    ACCUSED   sum    in   temp 


.  <0Df2; 
»ENn.U 

.  <nnFA  } 

. mupePj 

*n3: 
Tr^Pl 5 

J  P.  i 

Tfmpi ; 


'-en    FlMTS 


T  >  E  T  A  M  (  2  )     ,  ♦  1  } 
jump  LOOp.iS 


LOLC?) 

L0EE1 
Ll  HOl 

LbA 
SbpN 

STA 

LLLf  3"! 

2LPTA(3) 

L^  A 

StPN 

?TA 

L>A 

:  Lr  lc 2 •» 

LLEE1 

C**LL 

LIS 

L  I  T  (  3  } 

LP  A 
DVRN 

STA 

STA 

SLITC?) 

LUAD(?) 

LILC2) 

CAND(?i 

LDEE1 

LUA 

MLRN 

STA 

LUC?) 

CL*PC(?) 

LDEE1 

CUFAl 

STA 


»  <r>nrl  J 
so? » 

.  S  A  V  F  0  t 

0  (  0  )  ! 

T  F  mp  j 
0(0)  J 
•    iflpF  2  J 
,  \j  F  W  o  0  W 

T  r  M  p  i  5 
1  ( i'  •> ; 

Ii(l)U 
.  n>EP  ; 

s  a  R  T  6  ft  ( 
s  a  ; 

■  l  .o  j 

SC3| 

*$! 

0(0)) 
Temp  > 
«  T  E  m  D  * 
SC.3: 

, 'iHToEO 

unnFll 

sc2: 

0(0)| 

0(0)| 
. nODFI  I 


x  MnnE    pattern    rnR    LAST 

*  JUMP  If"  .K!V)E2  p-ALSE 
%  AMP  WITH  MQDEA 

%  AND  WITH  MODES 

*  NFW  MpDE  PATTEPN 
I  SUM  RESULT  IN  f,R 
5  WTTH  TEMPI 

X  SAVE  IN  TEMPI 


*  LOOP  K^E  »;HEM  J  =  fC?*I"l 

%  CONTINUE  LOOP  wHFN  [\"%^2<l"\ 

J  RFSET  MppF  PATTERN 

95  S  F  S  E  T  A  C  A  P  0 

*  SUBTRACT  SUM  (JN  TEMP) 
%  FPOh  I-TH  RPw  ANn 

?  STORE  hAC*  IN  I-TH  ROW 

i  MnPE  PATTERN  FPP  PARTIAL  ROW 

t  %  Jll*-  P  IF  MPPE2  FALSE  ' 

%  S '  I P  R  T  R  A  c  T  TEMPI 

*  FPPP  F£ST  OF  I-TH  ROW 

if.  STORE  MCK  IM  I-TH  ROW 

*  F  T  P  f,  T  PART  OF  ROW  TO  R  "•  A 

*  PATTEr-N  F0q  PEO 

*  f J m L  Y  PLO  ON 

W*  FIND  SORT  OF  FIRST  ELEMENT 


*  ro  SAvE  time  latEp  ppE«TE  recipro^al 
%  Gr  FOOT  (rAN  multiply  3ATmtp  THAN 

*  otvide  ^hfn  forming  d 

*  store  reciprocal  of  po"t  in  row  i 


%    MnVF  RESULT  OF  SORT  TO  SC3«1 /LC I • I  1 


Sn2, 
0(0)  J 


%  AmD  WITH  CURRENT  PATTERN 

*  PE-S  1  TO  MlN(f3#N-D 

4  LOAD    RCiA    WiTM    fOLD    ROk     I    - 

%  MnLT    Pt     1/Lt I'll 

%  REPLACE  A  BY  NEW  ROW  OF  I 


%    REVERSE  MPpE  PATTERN 

*  Zepo-fill  REST  OF  ROW 


SUM) 


LLL(?)    .MQDE2I 

ZLRTA(?)  »CKaC1;    X  JUMP  IE  MODE?cO  (I  CEQ  N-ffl) 


00028 
00028 
00028 
0002H 
00029 
00029 
00029 
00029 
00029 
00029 

nn029 
no  0  29 
00029 
00029 
00030 
000  3  0 
0  0  0  3  0 
00030 
0O030 
00030 
00030 
10030 

on03n 

00030 
00031 
00031 
0O031 
00031 
0O031 
('0031 
0  0  0  31 

^poii 

00031 
P0031 
00032 
00032 
00032 

0O032 
00032 
00032 
00032 
00032 
00032 
00032 
0OO33 
00033 
00033 
00033 
00033 
00033 
00033 
00033 

00033 
00033 
00034 

0003a 

0003« 


^no???? 

700???? 
«00???? 
900???? 

ooo???? 

100???? 

?oo???? 

300???? 
U  0  0  ?  7  ?  ? 

500???? 
600???? 
700???? 
800???? 
900???? 

o  n  o  ?  ?  ?  ? 

1 00???? 
200???? 
300???? 
ft  00???? 

soo???? 

60!)???? 
70"???? 
BOO???? 
9OO077? 

OOP???? 
1 00???? 
200???? 

30P???? 

'400???? 

5  00???? 
600  ???? 
700???? 

son???? 
poo????- 

OOP???? 
100???? 
?0U??7? 
300???? 

n  o  o  ?  ?  ?  ? 

500???? 
600???? 
700???? 
8  00???? 
901'???? 
000???? 
1 00???? 
200???? 
3  0  0???? 

aoo???? 

500???? 
600???? 
70U???? 

*no???? 

900???? 
000???? 
100???? 
200???? 


k2 


CKACl  » 


LI 
LI 

ML 
ST 
CO 
LU 
f.L 

si 

TX 
JU 


A  IS  ^ 
THUS  * 

1. 

?. 


% 
* 
X 

I 
X 
X 

SOLVElTjtF 

C 

s 
c 
s 

L 
S 
L 

L 
L 
L 

L 
l» 
L 
S 
.  K 

X 

L 
C 

X 

forlouf t »c 

L 
X 

L 
C 

t 
L 
S 
S 


EEl 
A 

RN 

A 

MPC(?) 

LEI 

FA  | 
A 

LFMM  ) 
MP 

OW  nEC 
i     HAvt 

soivr 
solv?: 

it  li 

ETE  E". 
LCC  1  W 
TLCl  l 

LRA» 
TA  TtW 
CL(  1  1 
TL(J  ) 
cEEl 
OL(?1 
OA 

Dl(C) 
P  A  0  (  0  ) 
LRN 

cccn 

TOpEfU 
TL 

DL(1) 

ADO( 1 ) 

ADD(O) 
FSSf  A( 

DL(3) 

SHR( 3) 
OR(3) 
PEE1 
T  L  C  3  1 
TR 


*C2| 

1(0)1 

1(0)  J 
I 
SC2» 


*  PjCK  I'P  REST  or  HOW 

*  MhLT  hy  1/LCI»I1 

*  STORE  bACK  IN  ROW  I 

X  REVERSE  MODE  PATTERN 


1(0),  %    ZEPC-FlLL    REST    OF    ROW 

N<Tpn>J       *Ju^'F    IF     I*SC1<N-1 

DMpn?FD  TMTn  I  AKD  LT(L  TRANSPOSE* 

I  *  LT  =  F.  LFT  (LT)X»y  ANO  THEN 

LY-F  rOR  Y  (FORMRO  SUBSTITUTION) 

(LT)X*Y  tor  x  (PAcK     SuesT I T  uT I  ON ) 


PR. "Ft  SETE1  E.^D.El 

,MPnE2»  %  sM  hOpE2  TO  FALSE 

P;  STA  TrKPl| 

,PFO;  *  MPPE  PATTERN  FOR  PEO 

.MpnElj  *  S/VE  TC  INDICATE  CURRENT  f^ODF 

Sfl  J  ?  Tl.RN  UN  PEO  OKLY 

.ftppfc;  *  atpr    df    ^composed    MTV 

0  (  ?  )  '  *  =  1  /  L  [  0  »  0  ] 

.AP&FI  X  j,(.(.sADUR    OF    F    ARPAv    (*X«F) 

*C1  I  T  %f  ) =Fl C  1 

*  c  j  i  *  YrojsFtovifo,oi 

*Aj  *  "CVE    TP    ACAR    1 

)     SCl>  »  nvERir-^lTF    F[Oj     wITM    YTO] 

$A.1i  *  Rf-UTF    TO    PE     1     -    STORE    IN    TFMD    LAjER 


,svi  onp;x-  L^r.p  counter 

.[IMF;    »  LOUP  FPOM  1  TO  N-l  (.flsROW  INDEX) 

.HKFI    X  PUMP  ACA»  0  in  GET  NEVT  El FMEnT  pF  F 
1)  SD2»F-nPMULT;  %     IF  I  GTP  63  JUM* 

.MpnElf  X  I  LEP  63 

U  *  GET  CURRENT  MpDE  PATTERN 

.hpnEi  j  >  Af'c  A[>U  T-Th  P.E 

?C3S  f  PE-S  0  THRU  I  ON 

,MnnEl»  '*  SAVE  PATTERN 

TEMP;  *  HJTH  PE  1  ON  f,R  CAN  U*ti     RF  STOREO 

V  Tf  '»  Pe(0,Y[  i-n,YtT-?3...vro)»o,..o) 


FORHULT, 


CLC(?)J 

CAOO(?) 

LOA 

SBM 

ApM 

STA 

LCLr?^ 
LDEE1 


sol  t 
«C?I 

PPNl 
«Aj 

RG*1  J 


X  RTA=I 

%  RGA=T-PFNs(I»!-l,...0) 

%  nOuPLE  VALUE  sINrF  2  *0ws  OF  PEM 

«  APf  PEUIJTREO  PER  ROW  ^F  A/L 

«  S/VE  FUR  1ST  PGX  PATTERN 


.wpoElt  ?  TURN  ON  PE-S  ACCHPO  Tn  MflOFl  PATRN 
*r?i     X  IF  I>63  ALL  PE-S  ON 


0003*300???? 
0p03<UOO?7?? 
COO345OO777? 
00034600???? 
00034700???? 

oco34boo???? 

00034900???? 
00035000777? 
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